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laminar  flow  regime.  The  time  dependent  components  at  the  forcing  frequency 
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momentum  equation  in  a  time  dependent  flow  requires  a  force  triangle  to  be 
maintained  at  any  instant.  The  triad  of  forces  are:  pressure,  inertia  and 
shear.  All  terms  of  the  force  balance  equation  were  measured  independently 
providing  a  good  check  of  data.  The  measured  turbulent  characteristics  of 
the  flow,  including  the  RMS  values  of  the  velocity  fluctuations,  Reynolds 
stress  and  short-time  power  spectra  are  dependent  on  the  phase  of  the  forced 
oscillations.  The  radial  distribution  of  HjLe  phase  angle  of  velocity  is 
qualitatively  different  in  laminar  and  turbulent  flows.  In  order  to  explain 
this  difference  the  concept  of  a  relaxation  time  of  the  turbulent  flow  was 
employed.  Transitional  flow  in  a  boundary  layer  is  also  briefly  discussed. 
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The  project, 
fully  turbulent 
complex  process. 

Research  is  carried  out  simultaneously  in  a  wind  tunnel  where  tran¬ 
sition  of  a  laminar  boundary  layer  is  extensively  investigated,  in  a 
pipe,  and  most,  recently,  a  large  two  dimensional  channel  facility  was 
added.  The  ability  to  compare  data  obtained  in  different  geometries, 
Reynolds  numbers,  and  pressure  gradients  enables  us  to  take  advantage  of 
the  unique  properties  of  each  facility  and  hiv<  an  overview  of  the  univ¬ 
ersal  aspects  of  the  problem. 

For  example:  transition  in  a  boundary  layer  takes  place  in  a  form 
of  "spot"  bounded  on  one  side  only,  by  a  solid  surface.  The  spots,  thus 
grow  spatially  in  all  three  directions  forming  complicated  three  di¬ 
mensional  structures.  The  pressure  gradient  over  the  surface  is  easily 
controlled  and  may  be  made  to  vanish  on  a  flat  plate  at  no  incidence, 
yet.  ‘here  is  no  control  over  the  growth  of  the  Reynolds  number  and  the 
local  thickness  of  the  boundary  layer.  In  a  pipe,  on  the  other  hand,  a 
transitional  turbulent  structure  "a  puff"  is  limited  by  the  diameter  of 
the  pipe  and  can  only  grow  in  the  streamwisc  direction,  the  Reynolds 
number  is  easily  controlled  by  controlling  the  mass  flux  but  the  pres¬ 
sure  gradient  in  a  given  pipe  at  a  fixed  Rc  depends  on  the  s  ructure  of 
the  flow,  and  is  thus  a  dependent  variable.  Although  the 


entitled  "On  the  relationship  between  transitional  and 
shear  flow"  is  concerned  with  various  aspects  of  this 
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•jtiscffibl  e-averageu  ‘uroulent  struct  ure  in  transitional  pipe  flow  is  axi- 
■'•ym;ne‘  rlc  ,  then-  is  lil.tie  evidence  that  this  is  a  dominant  mode  of  the 
transition  process.  The  two  dimensional  channel  facility  was  construct¬ 
'd  in.  cause  1'  provides  a  link  between  1  lie  boundary  layer  and  the  pipe 
*  low.  Al*  hough  'll-  pressure  gradient  In  a  channel  in  a  dependent  vari¬ 
able  the  transitional  structure  may  take  a  form  of  a  boundary-layer-spot 
if  generated  locally  by  a  point- source-di sturbance ,  or  a  "puff"  (or 
"slug")  if  generated  by  a  two  dimensional  disturbance  spanning  the  chan¬ 
nel.  However,  even  three  dimensional  disturbances  set  in  channel -flow 
can  only  propagate  in  two  directions  thus  theoretical  and  experimental 
analysis  of  the  proolem  is  simpler  then  in  the  boundary  layer. 

The  resulting  turbulent  flow  in  pipes,  channels  and  boundary  layers 
is  of  course  quite  similar.  The  viscous  sublayer,  the  logarithmic 
law-of-tne-wall ,  the  mean  turbulent  intensities  and  even  the  wall- 
streaks  and  large  coherent  structures  are  ail  alike.  This  implies  that 
only  those  features  which  are  common  to  all  geometries  are  of  paramount' 
importance  in  the  final  transition  process.  The  overall  research  pro¬ 
gram  on  transition  at  the  School  of  Engineering  in  Tel-  Aviv  was  guided 
by  this  view.  During  the  CY  lydo-tydl  w  >rk  was  done  on  the  Internal 
structure  of  the  boundary  layer  spot  (some  of  the  experiments  were  car¬ 
ried  out  at  the  University  of  Southern  California);  on  the  interaction 
of  spot 3  and  on  the  effects  of  favorable  pressure  gradient  on  spots.  In 
pipe  flow,  the  axl symmetric  mode  of  transition  was  investigated  by  mak¬ 
ing  d  simultaneous  measurements  at,  different,  av.1mut.hal  locations,  but 
the  bulk  of  tht3  years  scientific  report  i3  concerned  with  the  pulsating 
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pipe  flow.  The  introduction  of  harmonic  pulsations  to  pipe  flow  was 
proposed  by  us  in  1 y 72  but  was  not  executed  because  of  experimental  dif¬ 
ficulties  and  a  lack  of  an  adequate  data  acquisition  system.  The  pro- 
fa;;!  was  reestablished  in  1973  after  the  relationship  between  puffs  and 
iui'n  was  understood  (iiubin,  Haritonidis  and  Wygn. inski  ,  1930).  In  the 
a' inched  thesis  by  L.  Sheraer ,  pulsating  laminar  and  turbulent  pipe 
flows  are  discussed  in  detail.  The  response  of  the  system  which  con- 
1 ains  a  plenum  chamber,  a  piston  pump  and  a  pipe  to  controlled  slnosoi- 
d.tl  oscillations  was  investigated,  as  was  the  effect  of  the  finite  velo¬ 
city  of  sound  on  the  propagation  of  disturbances  in  a  long  pipe.  There 
are  marked  differences  between  the  structure  of  turbulence  in  a  pulsat¬ 
ing  pipe  flow  and  in  "steady"  fully  developed  turbulent  flow,  these 
differences  are  associated  with  the  memory  of  y-ig  turbuleqV  structure 
and  point  again  to  the  need  of  considering  relaxation  time  in  modeling 
turbulent  flows.  A  portion  of  the  thesis,  concerned  with  the  differ¬ 
ences  between  laminar  and  turbulent  pulsating  flows  and  the  modeling  of 
relaxation  times  using  complex  variables  to  denote  phase  relationships, 
will  be  presented  at  the  1981  Davis  conference  on  turbulent  shear  flows. 
The  effects  of  harmonic  pulsations  on  relam i nar I /ad ton  of  turbulent  flow 
and  on  transition  to  turbulence  will  be  considered  during  the  coming 


year  Cl  81-82. 


ON  THE  SUBSTRUCTURE  IN  A  TRANSITIONAL  SPOT 


By  1.  Wygnariskl 


Department  of  Fluid  Mechanics  and  Heat  Transfer 
School  of  Engineering 
Tel-Aviv  University. 


1)  INTRODUCTION 


Photographs  of  transitional  boundary  layer  spots  (Gad  el  Hak 
Blackwelder  and  Riley  (1980),  Cantwell,  Coles  and  Dimotakia  ( 1 97^ )  and 
Matsui  (1980))  reveal  that  the  spot  contains  numerous  eddies.  In  fact 
at  sufficiently  high  Reynolds  number,  the  interior  of  the  3pot  is  indis¬ 
tinguishable  visually  from  the  fully  developed  turbulent  boundary  layer. 
Longitudinal  streaks  aligned  roughly  with  the  direction  of  streaming 
figure  prominently  in  every  plan-view  photograph  of  the  spot  and  the 
boundary  layer.  Yet,  the  subtructure  within  the  spot  can  not  occur  en¬ 
tirely  at  random  in  view  of  the  universality  of  the  shape  of  the  spot, 
and  its  linear  growth  in  the  3panwise  and  streamwise  directions. 

Ensemble  avaraged  data,  conditioned  on  the  pertulbation  generating 
the  spot,  does  not  reveal  a  dominant  structure  within  the  spot  itself 
(Wygnanski,  Sokolov  and  Friedman  (1976),  Cantwell,  Coles  and  Dimotakis 
(1978)).  Consecutive  spots  may  uiffer  slightly  from  one  to  another  in 
their  shape  and  their  celerity  which  results  in  desynchronization  of  the 
acquisition  process  and  masks  the  internal  structure  deduced  fr'vs  1  he 
averaged  data.  Thus  the  velocity  field  reconstructed  from  ensemble 
averaged  data  would  lead  one  to  believe  that  the  spot  consists  of  a  sin¬ 
gle  large  coherent  eddy.  The  apparent  eddy,  however,  does  no*  scale 
correctly  with  any  characteristic  boundary  layer  1  e.ngi  h- scale  because 
the  spot  can  becomes  as  large  as  the  facility  In  which  it  is  generated. 


I'age  i 

Although  the  overall  dynamics  of  the  spot  is  of  some  practical  interest 
as  it  may  help  in  predicting  skin  friction,  heat  transfer  and  noise  gen¬ 
erated  in  the  transition  region,  it  can  hardly  explain  the  dynamics  of 
the  transition  process  itself. 

It  was  previously  observed  (Wygnanski,  llarl  t onl  d i  3  and  Kaplan  1  y  f 9  i 
that  a  wave  packet  trails  the  spot.  The  precise  relationship  between 
the  packet  and  the  spot  is  not  fully  understood,  particularly  at  low  Re¬ 
ynolds  numbers,  but  on  numerous  occasions  the  wave  packet  Dr’oke  down 
generating  a  new  transitional  spot  which,  at  least  initially,  did  not 
manage  to  catch-up  with  the  parent  structure.  The  breakdown  of  a  wave 
packet  into  structures  resembling  hair-pin  eddies  suggest s  that,  the  spot 
may  contain  a  fairly  orderly  array  of  such  eddies.  incipient  spot s  In 
favorable  pressure  gradient  (Wygnanski  !9d0)  contained  initially  three, 
then  five  distinctive  eddies  arranged  in  a  /\  formation.  The  number  of 
these  eddies  increased  for  a  while  with  down-  t ream  distance,  but  later 
they  merged  to  become  indistinguishable  further  downstreams  as  a  result, 
of  the  averaging  process. 

Thus  in  order  to  map  the  structure  of  a  "typical"  spot,  more  sophis¬ 
ticated  averaging  techniques  have  to  be  employed.  In  the  following  dis¬ 
cussion  few  attempts  are  made  to  pry  the  information  from  slreamwi se 
velocity  measurements  made  with  a  normal  hot-wire  rake  whien  supplies 
instanteneous  velocity  information  across  the  entire  boundary  layer. 

Measurements  were  made  in  the  absance  of  pressure  gradient  at  a 
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free  stream  velocity  of  5m/ sec  and  a  perturbation  located  at.  He  =  750 

o 

based  on  the  local  displacement  thickness  of  the  laminar  boundary  layer. 
The  experimental  facility,  the  Instrumentation  arid  the  data  acquisition 
system  were  discussed  previously  (Wygnanski,  Haritonldis  and  Kaplan 
’979)  and  will  not  be  repeated  in  the  presold  -on* ext . 


2)  ON  THE 


lii  A  3P0I 


An  example  of  !0  simultaneous,  at ream w i se ,  velocity  perturbation 
signals  ,  observed  during  a  passage  of  a  spot  is  shown  in  fig.1.  The 
abscissa  in  this  figure  is  time  while  the  ordinate  corresponds  to  velo¬ 
city.  The  numbers  opposite  each  trace  correspond  the  dimensionless  dis¬ 
tance  of  ci  given  wire  from  the  surface  of  the  plate,  expressed  in  terms 
of  tne  local  laminar  boundary  layer  thickness,  <T«-  .  One  may  clearly 
discern  large  velocity  fluctuations  which  arc  coherent  ocross  most  of 
idle  laminar  boundary  layer.  A  characteristic  frequency  associated  with 
these  fluctuations  corresponds  to  Uie  most  amplified 

"Tollmien-Schl ichting  frequency  appropriate  'o  inis  case  ( aec  also  Wyg- 
nanski,  Haritonldis  and  Kaplan  1979).  Power  spectra  averaged  over  200 
events  confirm  the  existance  of  energetic  fluctuations  at  tne 
Tol imien-Schl ichting  frequency.  The  spectral  peak,  however,  is  not  very 
strong  relative  to  the  background  turbulence,  and  it  is  not  clear  at 
this  point,  how  it  may  be  used  In  reconstructing  the  flow  field  associ¬ 
ated  with  a  typical  (or  a  most  probable)  spot  . 

After  examining  numerous  velocity  records  a  simple  criterion  was 
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de*  iced  which  classifies  the  record  according  to  easily  dl  st i ngui  siiabi  e 
number  of  large  coherent  fluctuations.  The  velocity  records  are  first 
low-pass-filtered  digitally  in  order  to  remove  from  each  record  the  ef¬ 
fects  of  the  small  scale  turbulent  fluctuations,  'The  filtered  record  is 
■.vanned  and  the  Mines  at  which  the  velocity  drops  to  a  minimum  locally, 
are  recorded  and  stored  in  memory  as  indicated  by  the  vertical  lines  in 
figure  2.  The  local  minimum  velocity  has  to  pass,  an  arbitrarily  deter- 
mind  t.reshold  level  relative  to  the  neighboring  maxima  before  being  re¬ 
corded.  The  Ireshold  criterion  eliminates  small  amplitude  fluctuations 
from  the  count  but  introduces  an  element  of  subjectivity  to  the  prneb- 
d  urn .  it  was  empirically  established  that  Wu:  following  il  soucsion.  {  s 
•i. sensitive  to  the  choice  of  the  t reahoid  level  .  For  example  hal fing 
'he  trcshold  level  would  result  in  accounting  for  one  additional  minimum 
it  rather  than  5)  in  the  upper  velocity  li  hoe  down  in  fig.  2. 

Adopting  this  procedure  for  each  measuring  station  enables  the 
classification  of  spots  according  to  the  number  of  large  eddies  which 
they  contain.  Some  typical  probability  density  distributions  showing 
'li'  most  probable  number  of  e-Jdlcs  >>u  the  plane  of  symmetry  of  the  spot, 
at  a  given  distance  from  the  perturbation  are  shown  in  fig.  j.  The 
most  important  conclusion  from  these  histograms  is  that  '.he  average 
number  of  large  eddies  in  a  spot  ]■  small. 

For  example;  measurements  m.we  375mm  downstream  of  the  perturba¬ 
tion  (975mm  from  the  leading  edge  of  the  plate)  indicate  that  the  aver¬ 
age  number  of  large  eddies  in  a  spot  is  4.0  and  the  most  probable  spot 
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detected  on  the  plane  of  symmetry  between  ()..:$  y/^  ^  0.8  ■;<•  n*  uin-d  v 
large  eddies.  The  most,  probable  number  of  eddies  n.-ar  '.he  surface  ana 
at,  y/g  >0.8  13  slightly  lower  (a!  y/ju s'  1  '.he  in-  at.  probable  number  1 :: 

h  >  because  the  Initial  breakdown  to  turbulence-  o.-.cirs  at  O.j *  y/j;t  £*  0.8 
( ,:c«?  ai  £>•»  Kovaaznay  Konjoda  and  Vasudova ,  '  yt>_: )  mi  some  of  tne  eddies 
tie.  not  arrive  perhaps  at.  the  inaesuring  station. 

The  eddy  detection  scheme  shown  in  fig.  8  renders  not  only  the 
number  of  eddies  in  a  given  realization  but  also  the  time  of  their  arri¬ 
val  at  tiie  location  of  measurement .  These  time 3  wore  stored  and  tne 

piobabilily  density  oj  their  dir.trlbul  ion  was  i:i;lne ■)  .  it.  appears  t.na' 

t.no  eddies  arrive  at.  preferred  times  at,  the  measuring  station  (fig.  41 . 
Thus  if  one  select.s  only  those  spots  wni.h  e<-tit  a  in  a  given  number  of  ed¬ 
dies  tne  histogram  showing  their  :-m-  of  urrivtl  <ii  1 1  contain  i  icntloai 
number  of  peaks;  each  peak  c>>rr-'  so-  nd  a  to  the  pt'el  erent,  i  al  time  of  ar¬ 
rival  of  the  specific  eddy.  Sine-,  idle  probibility  density  distribution 
of  the  number  of  eddies  in  a  spot  is  narrow,  by  choosi ng  only  those 
spots  which  contain  the  most  probable  number  of  eddies  for  further  ana¬ 
lysis  are  selects  between  25%  -  ’>>■'%  of  the  t-  '.al  number  of  events;  such 
i  large  fraction  justifies  the  subdivision  of  transitional  spots  into 
groups  containing  an  equal  number  of  eddies. 

The  abvementioned  information  can  be  used  obtain  the  detailed  flow 
field  in  a  spot  which  contains  a  proscribed  number  of  eddies.  The  pro¬ 
cedure  is  illustrated  schematically  in  fig.  b  showing  if:;  inherent,  ad¬ 

vantages  in  mapping  the  flow  fiele  d  a  "young"  spot  containing  j  eddies 
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onl  y : 

(i)  The  number  of  eddies  in  a  spot  is  detected  and  i  simple 

:  n  sumb  1  e- a  v  c raged  veloei'.  y  condition-  I  -u  ’  |  vi  1  -u  n  it .  -  >ti  i  ■:  :  >j 

euiated  (fig.  5a).  .Since  only  the  :;p -1  •-  imhi  .lining  i  w-.-r  : 

selected  for  the  averaging,  the  ensemble  averaged  velocity  contains 
i  local  minima.  it.  is  clear,  however,  that  the  amplitude  of  tin.  so 
minima  is  not  representative  of  the  velocity  record  in  a  single  re¬ 
alization  drawn  io  the  same  scale  in  fig.  5 f . 

(ii)  Each  realization  is  shifted  in  tiino  in  order  to  align  i  he 
location  of  its  last  minimum  before  averaging.  lho  ensemble  aver¬ 
aged  velocity  resulting  f ran  the  alignment  procedure  Is  shown  in 
fig.  5b.  The  amplitude  of  the  trailing  fluctuation  in  '.he  ensem¬ 
ble  is  enhanced  by  the  alignmont.  process  .m  l  is  compurcahlc  to  (  n<‘ 
amplitude  of  the  single  realization  shown  in  fig.  51.  1'he  neigh¬ 
bouring  minimum  was  hardly  affected  by  the  alignment  wliile  the 
first  minimum  disappeared. 

( ill  Repeating  the  alignment  procedure  for  the  first  and  i n< 
second  minimum  enhanced  '.he  amplitude  of  the  particular  minimum 
concerned  (fig.  5c,  5d) . 

(iv)  A  composite,  ensemble-averaged  velocity  perturbation  re¬ 
cord  is  generated  by  subdlvidiig  the  overall  temporal  record  into  4 
sections  (fig. 5c),  This  subdivision  sugg.  d  :  ini'  .a.-h  l  u-g.-  eddy 
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has  a  limited  zone  of  influence.  The  velocity  per'  urtiai  i  on  in  re¬ 
gions  (!)  and  (4)  a re  under  the  total  influence  of  minima  A  and  C 
respectively.  The  velocity  perturhat  ion  in  regions  (2)  and  (3) 
have  heeti  weighted  proporf  i  n  ■  1  1  y  to  in-'  l  i  • '  mee  •  >  1  ■  p  n  I) 

from  points  A,  B  and  C  respective!  V  :>  Ion-  averaging.  Thus  ‘re 
velocity  preturbat Ion  at  points  A;  H;  C;  are  identical  to  the 
velocity-perturbations  snown  in  figs.  5d,  be,  ub  respectively. 
Tiie  velocity  in  region  2  is  a  weighted  average  of  tiie  records  shown 
in  fig.  5c  and  5d  respect.1  vel y  while  the  velocity  in  region  3  is 
the  weighted  average  of  records  5c  and  51). 

This  procedure  can  be  applied  to  any  number  ->f  eddies  -xis'itig  in  a 
most  probable  spot  at  a  given  coordinate  in  spice.  The  ultimate  purpose 
of  the  process  is  to  reconstruct  a  detailed  j  dimensional  flow  field  in 
*  ne  most,  probable  realization  from  4  large  number  of  point  measurements, 
icoeptiiig  the  fact  that,  consecutive  realizations  are  only  broadly  simi¬ 
lar  .  The  success  of  the  procedure  is  determined  by  comparing  the  recon¬ 
structed  velocity  perturbation  contours  in  a  op"'  with  'he  contours 
measured  by  an  array  of  wires  during  a  passage  of  a  single  realization. 
The  comparison  is  made  al  the  largest  distance  from  t ne  location  of  the 
perturbation  for  which  data  is  available  (i.e.  at  X-Xpert  =  375mm).  At 
this  distance  the  spot  is  already  quite  large  and  may  be  considered 
fully  turbulent  (Rex^£  >  3 • 3* ! 05 or  He  S*  ~  "**  ).  The  comparison  thus 
repr3ents  the  most  severe  test  of  the  procedure  based  on  the  available 
da'  a  . 
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A  simply  ensemble-averaged  velocity  perturbation  record  is  shown  in 
fig.  6a.  The  velocity  perturbation  is  entirely  positive  near  the  sur¬ 
face  and  entirely  negative  at  large  ii  stances  from  1  lie  wall.  l'his  re¬ 
cord  ,  however,  differs  somewhat  from  similar  dot  »  reported  by  Wygnanski 
Haritonidis  and  Zilberman  (1980)  because  trie  me  inurements  were  made  at 
much  lower  Re.  The  most  obvious  difference  concerns  the  existence  of 
two  minima  at  an  intermiate  distance  from  the  wall,  (0.25<y/j  <0.b) 
which  are  marked  by  arrows  in  fig.  ba.  Another  difference  is  concerned 
with  'he  existence  of  a  positive  velocity  perturbation  near  the  leading 
edge  of  '.he  spot  which  extends  io  y/gu  ~0.b,  while  it  no  longer’  exists  at 
y/fL  s0.2*>  at  higher  He  (see  Wygnair  k:  Haritonidis  and  Zilberman  1 980 ) . 

The  same  data  was  reprocessed  according  to  the  procedure  outlined 
earlier  and  plotted  in  fig.  6b.  This  data  represents  a  velocity  per¬ 
turbation  resulting  from  the  sole  uiflueuce  of  tin  large  eddies,  in  the 

most  probbable  realization,  on  the  plane  of  symmetry  of  the  spot  at  a 

given  Re,  and  X.  The  most  probabl  ■  spot  contains  5  large  eddies  at 
0.  15  $  y/Su  $  0.6;  4  eddies  near  the  surface  and  a1  0.65  <.y/g,S.  0.85  and 

only  j  eddies  at  y/j^l.  The  perturbations  sh"wu  in  fig.  bb  contain 
detailed  Information  ab.iut  the  interior  structure  of  the  flow  which  is 
not  visible  in  fig.  6a.  Heavy  filtering  of  t  tie  la' a  in  fig.  6b  repro¬ 
duces  the  regular  ensemble  averaged  vcloci'y  }*■  t  t  urbat  Iona.  Velocity 
perturbations  recorded  during  tin;  passage  of  '  u.  (irst  spot  In  the  ex¬ 
periment  are  shown  in  fig.  6c.  The  velocity  tan  ord  contains  high  fre¬ 

quencies  associated  with  small  turbulent  eddies  which  are  3imply  fil¬ 
tered  and  replotted  on  the  same  figure.  He-  filtered  signal  is  di s- 


Page  9 


placed  vertically  for  convenience  of  compart  son ,  and  although  It.  1;;  noi 
;  nsi*  ieal  to  the  most  probable  realization,  the  latter  appears  to  eon- 
•aln  the  most  relevant  information  which  is  missing  in  fig.  ha. 

Velocity  perturbation  eohnl  our:;  non  <:;po|ii  ;  i,f;  to  fig:;,  ha, bp;  are 
plotted  in  fig.  7.  The  vertical  scale  in  this  figure  corresponds  to  a 
distance  from  the  wall  while  the  horizontal  scale  is  time.  Using  a  re¬ 
presentative  convection  speed  it  appears  that  the  vertical  scale  Is  ap¬ 
proximately  stretched  by  a  factor  of  bu .  The  similarity  between  the 

velocity  perturbation  contours  recorded  during  a  single  event  and  the 
'•du-'od  most,  probable  perturbation  contours  shown  in  fig.  7h  is  impres- 

sivt,  while  the  contours  shown  in  figur-  contain  no  such  dot. ail. 

3)  IM  SPANWlSt:  STRUCTURE  M  A  TRANSITIONAL  oPyT 


There  are  objective  difficulties  in  obtaining,  reliable  data  near 
the  wall  using  a  rake  of  hot  wires  oriented  in  the  spanwise  (Z)  direc¬ 
tion,  because  any  small  surface  positioned  parallel  to  the  wall  may  pro¬ 
duce  lift  when  inserted  into  a  boundary  layer  resulting  in  the  shedding 
of  trailing  vortices.  One  may  not  overcome  this  difficulty  by  position¬ 
ing  the  rake  parallel  to  the  direction  of  the  mean  streamline  at  every 
level  from  the  surface  because  the  rake  may  still  lift  per i od 1  call y 
while  interacting  with  the  large  coherent  structures  in  the  boundary 
layer  and  shed,  in  addition  to  streamwi  ;ie  vorMeity,  ,i  starting  vortex 

which  is  equivalent  In  strength  to  the  bound  vortex  associated  with  the 

•  -* - *  • . -  — 1  ..iii.  ,  j.rli<e  of  wires 


1  i  ft . 


Thus,  all  present  experiments  were  made  with  a 
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which  are  normal  to  (tie  surface  ( see  figure  1  d  Wygria  .ski ,  Marl  toi;J  ji  a 
and  Kaplan  1979)-  The  flow  was  '  r  .ivirwl  from  the  plane  if  :iyinn"  tr  y 
outwards  at  Intervals  of  AZ=2ran  in  order  *o  provide  suffi.-ienf  spatial 
i  P  ion  of  1  hr  longitudinal  sin. s'  rue*  ure  . 

Measurements  in  favorable  pressure  gradient  (WygnarisKi  I960)  have 
shown  that  the  velocity  perturbation  at  a  given  distance  front  tne  wall 
is  not  homogenously  decaying  across  the  span  of  the  spot.  Positive  vel¬ 
ocity-perturbation  regions  near  the  sol  in  surface  w ere  separetou  t»> 
areas  of  decelerated  flow,  point. 4  ng  to  the  possible  existence  of  long*  - 
’  ud  itial  vortices.  The  number  of  t lie  stru.. •Lures  Increase  I  i u 1 1  '.a 1 1  y  in 

•  ne  !  wn  stream  direction  untill  t  n«  longitudinal  •  ubs'.ru  ji  ure  bee  ».jo  i.  r.- 

distinguishable  in  the  mean,  giving  an  appearen.o  of  a  single  largo 
w-rco- shoe  eddy.  it  was  then  realised  'ha*  ■  n-  individual  spots  vary 
considerably  in  width  ,  so  that  the  averaging  pr< m.-sa  includes  real i za- 
tions  which  do  not,  extend  as  far  in  7.  as  the  l  ine  of  measurement  and  uo 
n>  t  contribute  to  the  velocity  perturbation.  lin-  saw  data  was  then  re- 
pi  after  t  he  smaller  spots  w  re  exel  aU  d  t  s  - m  *  he  averaging,  and 

'he  ep, unwise  st  ructure  re-emerged  :n  the  mean.  tic.  pro  euur  is  biased 
'  .wards  the  larger  spots  but  it.  n:r;  'he  ad  van'  ire  I  sit'  i  ng  .  iu*.  in  •ele¬ 
ven.'  data. 

A  simple  ensemble-averaging  pr< ce  l  ure  eon  i  I '  '  •  ■  i > ■  d  t  o  hi.  perturb. i - 

*  Ion  was  used  in  t  he  present  experiment  1  o  e.|  ice  ■  ) i ■  ■  ui.inwi  ;»  st  rue!  ure 

1 ‘.(hun  downstream  of  the  origin.  The  incipim'  (fig.  dr  ;..nt  lin.- 

t  tin  r  longitudinal  structures  situated  in  a  A  Ioi'm'  fa,  j  ...  ,|  ,,  Wyg- 


Kigr i  I  I 

ri  <ns*i  19t)0).  The  spanwise  fil  stance  be' ween  t  ne  m  nters  of  these  s'.ruc- 
'  virea  near  the  surface  In  Verms  of  w.  1 1  coordinal  es  -9'>  (where 

'  >ie  friction  velocity  U*  was  estimated  from  fully  developed  turbulent 
boundary  layer  data).  The  number  <1..  s<#ii  is  ltd  •.  r*.  :?1  irig  ,  be.uuse  if  is 

>*fina*  i;ly  oqtnl  to  1  fie  •list  u.ee  <n<-;i:;ur--.-|  t>.  •*  worn  longitudinal 

freaks  In  a  fully  t.urbulent  boundary  layer.  Hen. siting  the  same  meas¬ 
urements  In  a  developed  turbulent  spot,  (at  x-xperj:  jOOmm)  nul  biasing 
'he  averaging  towards  the  larger  events  indicates  that.  the  spo'  contains 

Vi  structures  In  the  spanwise  direction.  There  is  an  increase  in  the 

number  of  the  structures  with  downstream  dl  s'  .anno  but  there  is  also  an 
’  a  'reuse  in  the  distance  between  ‘he  : enters  d  i  lja  .a.nt  structures.  It 
’.relieved  tnat  the  increase  in  the  number  •  .  ‘ met  urea  is  dominating 

'in  lateral  growth  of  t.he  spot. 
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1.  Velocity  signals  during  a  passage  of  i  lie  spot . 

.  .  i fi*;  del cet. i r>n  of  s'gni  ficani  minima  wi<  i.i \-  ■’  . 

I .  A  Histogram  showing  the  number  of  eddies  in  a  spot. 

A.  A  Histogram  showing  the  times  of  arrival  of  tno  large  eddies  ai  tile 

measuring  station. 

5.  The  eduction  of  the  large  coherent  eddies. 

(a)  Simple  ensemble  average.  ( b )  Ensemble  averaged  velocity 
aligned  on  the  most  rearward  e.iiy.  (>•)  Ensemble  averaged  velo¬ 
city  aligned  on  the  middle  eddy.  (  l)  Ensemble  averaged  veloci¬ 
ty  aligned  on  the  leading  eddy,  (e)  A  weighted  compos  te  aver¬ 
age.  (f)  The  velocity  in  a  single  realisa'  ion. 

o.  Velocity  records  in  a  spot. 

(a)  A  simple  ensemble  averaged  record.  (b)  A  composite  ensem¬ 
ble  averaged  record  showing  the  most  probable  distribution  of 
,  the  large  eddies,  (c)  The  velocity  perturbation  during  a  sin¬ 
gle  event. 

7.  Contours  of  velocity  perturbation  correspond  i  1 .  •  fig.b  a  ,  b, 
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The  ensemble  averaged  spanwise  structure  in  a  spH. . 

(a)  (X-Xpert  )=  150mm.  |  (b)  X-Xpert  =  300min .  (c)  A  schematic 

plan  view  of  a  spot. 
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ABSTRACT 


Pulsating  flow  of  air  in  a  straight  smooth  pipe  was  investigated 
experimentally.  Most  measurments  were  made  at  a  mean  Reynolds  number  of 
4000,  but  the  influence  of  Re  was  checked  for  2900<Re<7500.  The  period 
of  forcing  ranged  from  0.5  sec  to  5  sec  which  resulted  in  the  change  in 
the  non-dimensional  frequency  parameter  ot  =Rv  4.5  to  15.  Velocities 
at  the  exit  of  the  pipe  and  pressure  drop  along  the  pipe  were  measured 
simultaneously;  velocity  measurments  were  made  using  arrays  of  hot-wire 
anemometers  consisting  either  of  normal  wires  or  an  x  array.  Signals 
from  the  anemometers  or  pressure  transducers  were  digitized  and  pro¬ 
cessed  by  a  minicomputer  before  being  recorded  on  a  magnetic  tape. 

The  Introduction  of  periodic  surging  had  no  effect  on  the  time  mean 
quantities.  The  present  data  was  compared  in  detail  with  the  theory  of 
Uchlda  (1956)  in  the  laminar  flow  regime.  The  time  dependent  components 
at  the  forcing  frequency  were  represented  by  the  radial  distribution  of 
amplitude  and  phase.  An  integral  momentum  equation  in  a  time  dependent 
flow  requires  a  force  triangle  to  be  aaitained  at  any  instant.  The 
triad  of  forces  are:  pressure,  inertia  and  shear.  All  terms  of  the 
force  balance  equation  were  measured  independently  providing  a  good 
check  of  data.  The  measured  turbulent  characteristics  of  the  flow,  in¬ 
cluding  the  RMS  values  of  the  velocity  fluctuations,  Reynolds  stress  and 
short  time  power  spectra  are  dependent  on  the  phase  of  the  forced  oscil- 
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lations 


The  radial  distribution  of  the  phase  angle  of  velocity  is  qualita¬ 
tively  different  in  laminar  and  turbulent  flows.  In  order  to  explain 
this  dlfferenoe  the  concept  of  a  relaxation  time  of  the  turbulent  flow 
was  deployed.  A  staple  eddy  viscosity  aodel  for  time  dependent  flow, 
which  takes  into  account  the  "memory"  of  turbulence  is  proposed,  and 
numerical  solution  of  the  Navier-Stokes  equation  for  the  turbulent  pul¬ 
sating  pipe  flow  utilising  the  proposed  model  were  obtained. 
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Courier  transform  coefficients 
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= arbitrary  functions 
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slengtn  of  the  pipe 
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CHAPTER  1 


IMTROPUCTlOw  KhU  REVlUm  Of  LxTbnATUhb 


1.1  General  Introduction 


The  pulsating  pipe  flow  is  considred  a  "simple”  time  cicpcr.eer.t 
flow.  it  is  "simple"  both  spatially  and  temporally.  In  tne  aosencc  of 
swirl  the  pipe  flow  i3  two-dimensional  spatially,  anc.  in  tne  fully  -de¬ 
veloped  region  all  mean  quantities  are  functions  of  tne  radius  only. 
Nevertheless,  the  theoretical  solution  for  tne  proDiem  of  the  linear 
stability  of  the  stationary  pipe  flow  was  only  recently  obtained  (see 
Goldshtick  and  Shtern  (1977)  and  the  process  of  transition  from  tne  lam¬ 
inar  to  turbulent  flow  is  still  an  enigma  although  it  nas  oeen  investi¬ 
gated  experimentally  (kygnanski  and  Champagne  (1973),  Rubin,  wygnansai, 
and  Haritonidis  (1979)).  The  spatial  "simplicity"  of  the  flow  is  tnere- 
fo re  somewhat  delusive.  The  superposition  of  simple  harmonic  oscilla¬ 
tions  on  the  steady  mean  flow  adds  a  temporal  dependence  which  consider¬ 
ably  complicates  the  detailed  analysis. 

The  importance  of  studying  time  dependent  flows  in  general,  and 
pulsating  pipe  flow  in  particular,  is  obvious,  host  biological  flows 
are  pulsating,  may  be  beoause  the  peristaltic  pump  is  the  simplest  pump 
which  can  be  employed  by  a  biological  system.  Pulsating  pipe  flows  were 


i'age  j 


tnerefore  extensively  studied  oy  investigators  associated  nitn  medicine 
sr.o  life-sciences  (for  references  see  daro,  peaiey  ano  ocnroder  (15,70), 
nussain  ( 1 977 1 1  -  won-steady  flows  occur  also  in  many  engineering  appli¬ 
cations,  for  example:  tne  ciscnarge  of  any  piston  pump  is  pulsating, 
tnus  the  flow  in  an  intake  or  exnaust  manifold  of  internal  eomoustion 
engine  is  pulsating;  the  flow  in  hydraulic  lines  anu  control  systems 
often  pulsates  etc. 

1.2  ine  analysis  of  Periodic  uata 


‘me  periodic  nature  of  the  turbulent  flow  suggests  tne  decomposi¬ 
tion  cf  any  flovr  variable  gix.t)  into  3  components  (hussain  and  heynoids 
(1  Wvj) 


r (x,t)=g(x)+<g(x,  <p)>*g'  U,t) 


1.1 


vnere  six)  is  time  mean  value  of  tne  variable  g(x,t;  at  point  x 


_  T/2 

g(x)  a  liffl  1/1  [g(x,t)dt 
T-teo  J 
-T/2 


<s£ix,  <P)>  is  tne  contribution  or  the  periodic  part  at  a  pnase  angle  <J>  at 
tne  same  point  x  and  is  defined  oy 


<g(x,<p)> 


?  -  “ 
*  lim  1/ft  l  i&Ax,$  )-gix)) 

N-x»  i«l  1 
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ana  g’(x,t)  is  tne  random  nonperiodic  part.  experimentally,  the  value 
of  <g>  is  determined  oy  averaging  tne  aeta  at  a  fixed  pnase  angle  ana 
subtracting  from  the  result  the  time  mean  value  g.  The  flow  is  pulsat¬ 
ing  provided  the  quantity  <g>  does  not  vanish.  Prom  the  definition  of 
the  phase  mean  component  <g>  it  is  clear  that  the  following  relations 
must  hold  for  tne  temporal  mean  values 


g's<g>sg'<h>sO  1.2 

A  distinction  is  sometimes  matie  between  "pulsating"  anc.  "oscillating" 
flows;  the  former  term  implies  that  oscillations  are  superimposea  on  a 
non-vanishing  steady  velocity,  while  the  tern  "oscillating"  flow  refers 
to  g(>:)=0.  'ihe  phase  dependent  part  of  the  flew  <g>  is  referred  to  as 
tne  oscillating  part. 

assuming  tnat  the  flow  is  not  only  periodic,  out  also  Harmonic,  ana 
defining  the  phase  angle  $=  cot,  tne  periodic  component  may  be  represent¬ 
ed  oy  tne  reai  part  of  the  exponential 

<g(x,  $)>sne(gj(x}exp(i(  <ot+  tg>>  l.j 

wnere  g^(x)  is  tne  amplitude  of  pulsations  and  is  a  reference  pnase 
angle.  All  pnase  angles  referred  to  in  this  woric  are  measured  relative 
to  pressure.  Whenever  a  periodic  motion  is  not  simply  harmonic,  tne 
time  dependent  variable  can  always  be  expanded  in  courier  series,  so 
tnat  the  right  hano  side  of  tq.(lo)  becomes  tne  leading  time  dependent 


tern  in  the  expansion 


in  tne  analysis  of  the  narionicaliy  pulsating  flows  it  is  often  co- 
venient  to  use  vectorial  notation,  which  is  generally  aceepted  in 
electrical  engineering,  and  represents  graphically  the  periodic  compo¬ 
nent  of  the  variable  <g>  by  a  vector  of  the  iengtn  £  ana  angle  <p  re- 

f  & 

lative  to  a  reference  direction  (in  our  case  tne  pressure}.  This  pre¬ 
sentation  will  be  used  later  for  further  elucidation  of  pnase  and  ampli¬ 
tude  relationships  in  a  pulsating  pipe  flow. 


1.3  Governing  Equations 


In  the  absence  of  swirl  the  circumferencial  component  of  velocity  w 

vanishes  thus  only  the  axial  u  and  radial  v  velocity  components  has  oeen 
*1 

considered.  '  The  continuity  equation  for  tne  incompressible  flow  in 
cylindrical  coordinates  has  the  following  form: 

3u/3x+1/r  3/3r(vr)*Q  1.4 


Decomposing  the  axial  and  tne  radial  velocities  into  tnree  distinct 
components  according  to  (1.1),  and  averaging,  we  obtain  that  the  contin- 

i) 

The  full  set  of  equations  for  the  nonstationary  flow  in  the  cylindr¬ 
ical  coordinates  was  given  by  Cneng  (1971). 


uity  equation  in  tne  form  of  (1.4;  is  valid  for  each  of  tnrce  components 
separately:  time  mean,  periodic,  and  random. 


*ne  i.avier-stoKes  equations  for  axial  and  radial  velocity  compo¬ 
nents  (see  hinze  ( 1 V75 ) )  are 


3u  3u  3u 
37  +u  37  +  v  3?  = 


1  9£  +  V(iJi 

P  ox  \  2 


32u 

+  -  + 

2 


3r 


1  3ik 

7  3? 


3v 

3t 


♦u 


3v  1  3p  x 

v  37  =  '  p  37  + 


92v  1  3v 


1.5 


1.6 


decomposing  the  velocity  and  pressure  terms  in  (1.5;  and  (1.6)  ac¬ 
cording  to  Eq .(1.1)  and  averaging  with  respect  to  time  taxing  into  ac¬ 
count  the  definitions  (1.2)  yields: 


-  3u  -  3u 
uH  3r 


-  £  H  *  I  ST  '  T  ST  !•» 

OX 


-  9v  -  3v 

u37+v  37 


?ll*  <f7*7  5Fr5F)-FlFr(<v>I*v'i)-|r('<u><v>*u'v,)  l'8 


Equations  (1.7)  and  (l.d)  indicate  that  the  Reynolds  stress  in  pul¬ 
sating  turbulent  flow  can  be  regarded  as  a  sum  of  the  oscillating  ana 
random  parts.  The  Reynolds  equations  for  the  oscillating  component  may 
oe  deduced  by  phase  averaging  of  the  time  dependent  equations  (1.5)  and 
(l.b)  and  suostracting  from  the  result  equations  (1.7)  ana  (1.6),  res¬ 
pectively. 


men  tne  pipe  flow  is  fully  ee/eiopsu,  £.11  toe  time  ana  pnase  mean 
pare :.;eterr  .-.re  in^epenaent  ot  toe  atreamvisc  cocrc.in5.te  i;  tnus  it  fol¬ 
low?  free  tne  continuity  ec.u2.ti0n  t.'.at  tr.c  temporal  e.na  pnase  aver a*e>- 
r-2..iei  components  oi  velocity  vanish  in  the  fully  oeveiopea  region. 

Vs<V>»C 

The  Aeynoios  equations  for  tnc  stationary  ane  oscillating  components  of 
tne  axial  velocity  are  then: 

t's-1/p  3p/3x+1/r  3/3r(  vr3u/3r-ru’v'i  l.v 

9<u>/3t--1/p  3<p>/3x+1/r(  vr3<u>/3 r-r<u'v»>)  1.10 

ana  for  the  fuiiy  developed  laminar  flow,  equation  (l.jj  is  reuuccu  to 

3u/3t»-1/p  3p/3x+1/r  3/3r(  vr3u/3r)  1.11 

As  i  result  of  linearity  of  (1.11),  the  solution  is  a  superposition  of 
two  independent  parts:  tne  steady  parabolio  velocity  profile  calculated 
oy  roiseuille  and  the  time  dependent  osoillating  part. 

I.w  neview  of  Pertaining  Literature 


Velocity  measurements  in  oscillatory  pipe  flow  were  first  mace  oy 
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a isn:.rcison  < hi  t/Zii)  and  r>icnarason  and  Tyler  vHav).  Tney  found,  tnat 
in  pul  satin;;  laminar  flow  tne  aaxiaua  oi  tioe-me&n  velocity  occurs 
near  tne  wall  of  tne  pipe  ratncr  tnan  on  tne  centerline  as  it  does  in 
stationary  flow.  Tne  location  of  tne  maximum  velocity  approached  tne 
surface  witn  increasing  the  frequency  of  pulsation,  inis  phenomenon  was 
called  "annular  effect". 

1.4.1  Laminar  flow 

oexl  (153d)  gave  theoretical  explanation  to  tne  annular  effect  in 
laminar  flow  oy  solving  equation  (1.11).  The  lully  developed  laminar 
oscillating  velocity  profile  was  found  to  be  controlled  by  a  single  non* 
dimensional  frequency  parameter  a  *h/u>A> .  ror  low  values  of  tne  fre¬ 
quency  parameter  (  a«1)  the  velocity  profile  is  nearly  paraoolio  in  ac¬ 
cordance  with  the  instanteneous  value  of  tne  pressure  gradient,  for 
nigner  values  oi  tne  velocity  profile  deviates  from  the  Poiseuille 
form  and  the  annular  effect  appears. 

nooersley  (1555)  and  ucnida  (155o)  extended  tne  solution  of  sexl 
for  tne  velocity  distribution  and  calculated  the  phase  and  amplitude  re¬ 
lationship  among  the  pressure  gradient,  the  mean  bulk  velocity,  and 
shear  stress  in  oscillating  or  pulsating  pipe  flows.  The  paper  of  bohi- 
da  gives  the  solutions  in  tne  most  convenient  fora  and  this  trill  be  used 
as  the  main  reference.  Atabe*  and  Chang  (ISbl)  tried  to  solve  equations 
(1.5)  and  (l.o)  for  tne  developing  pulsating  flow  in  the  entrance  region 
of  the  pipe.  They  estimated  the  non-linear  inertia  terns  by  assuming 
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taea  to  us  pro,  crtio-u  1  to  tnc  iooc.i  rclus  oi  3u/3x  ant  to  t.ie  instan- 
tcneor.s  vain  a  t.'.c  velocity  entrance  at  t;.w  ul  1-:-e.  cos  a»t> ,  wnerc  ~  is 
t..e  value  ci  t.x  vr.ifcra  veiccit.,  at  tr.e  entrants  ox  t..c  pipe,  ana  a  is 

xv 

relative  xmpxitu...^  ox  velocity  pulsations,  me  results  are  not  present¬ 
ee  In  an  analytical  fora  out  a  proccuure  is  aevexepec  mica  provlues  a 
solution  for  any  given  set  ox  parameters. 

Linford  ana  hyan  Hvo£>l  measured  the  relation  oetv-cer,  tne  pressure 
gradient  ana  tne  flow  rate  in  an  oscillatory  iio  o*  a  mixture  oi  water 
uitn  glycol,  xiiey  used  only  flow  visualisation  tccnr.icues  for  measuring 
velocity  ana  tne  accuracy  ox’  tne  results  was  po or.  mey  concludes, 
novcver,  tnat  the  measured  velocities  agrees  with  tnc  theoretical  calcu¬ 
lations  of  womerslcy  within  tnc  estimated  experimental  error. 

tsnison  \  1 V- 70 >  and  Lenison  et  ai  U'.71v  repcrtcc  on  mcasurments 
made  in  pulsating  laminar  fxov  with  uirectionaiiy  sensitive  laser  vclo- 
cimcter.  iheir  apparatus  was  500  diameters  ion,;,  so  tneir  measurments 
coula  oe  mace  cotn  in  tne  entrance  region  ana  tnc  xuily  ceveiopea  region 
of  tne  pipe,  me  viscosity  of  tne  wording  fluid  was  v-.riaule  in  ora er 
to  proauce  tne  desired  flow  parameters  (mean  ney»iOj.as  nucuer  aoout  1000, 
anti  frequency  parameter  4£ft£o).  Xne  results  ootainca  in  tne  fully  ae- 
velopea  region  were  in  gooa  agreement  with  theory  out  in  tne  entrance 
region  aiscrepancies  were  founa  oetween  tne  measurements  ana  the  predic¬ 
tions  of  Ataoew  ana  uneng.  xnc  measured  velocity  profixe  cevclopeu  sig¬ 
nificantly  slower,  tnan  predicted  tneoretically. 


•*m*  WflfaMWi :  in. 


Omni  ct  al  U97b)  investigated  Doth  tneorctically  ana  experimental¬ 
ly  tne  laminar  pulsating  pipe  flow  in  air,  taxing  into  account,  tnat  tne 
fluia  is  3ligntly  compress! ole.  They  concluded  that  for  fl<~*«3a,  wnere 
a  is  tne  velocity  of  sound,  tne  radial  velocity  and  cross  sectional  var¬ 
iations  of  pressure  are  negligioly  small,  inis  result  implies  that  the 
compressibility  of  the  fluid  does  not  change  the  radial  distribution  of 
velocity  as  long  as  the  radius  of  tne  pipe  is  much  smaller  than  the 
sound  wave  lengtn  produced  by  the  oscillations. 

1.4.2  Turbulent  and  Transitional  flows 

One  of  the  interesting  features  in  tne  laminar  pulsating  pipe  flow 
is  that  it  nay  contain  temporarily  an  inflection  point  at  some  phases  of 
the  motion.  At  nigh  values  of  the  frequency  parameter*,  points  of  in¬ 
flection  exist  at  all  phase  angles.  It  is  well  Known  tnat  an  existence 
of  the  inflection  point  is  a  necessary  condition  for  tne  instability  of 
unviscid  flows  (ftaylelgh  criterion).  Sarpkaya  (19b6)  tried  to  correlate 
the  instability  of  oscillating  pipe  flow  with  tne  existence  of  an  in¬ 
flection  point  in  the  velocity  profile,  he  found,  that...  "the  flow 
has  maximum  stability,  when  the  duration  of  the  inflection  period  re¬ 
aches  about  53>  of  period  of  pulsations".  (Sarptcaya  (156b),  p.596).  It 

snoulo  be  noted  that  Sarptcaya  used  erroneously  a  plane  criterion 
2  2 

(d  u/dr  >0)for  instability,  i.e.,  instead  of  a  criterion  applicable  to 
cylindrical  coordinates,  namely:  d/dr(r  d/dr(ru))»0  (Batchelor  and  Gill 
(1962)).  Nevertheless,  these  errcr-  could  not  affect  his  conclusions 
significantly.  The  stability  of  bounded  periodic  flows  seems  to  have 
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little  in  common  witn  the  iD»ylei«n  criterion  for  an  ur.viscid  fluid,  in 
orntrast  to  free  snear  flows,  litre  mixing  layers,  wa*:es  and  jets,  wnose 
instability  is  governed  by  an  inviscid  process,  tne  appearance  of  an  in¬ 
flection  point  in  an  oscillatory  Dounoed  flow  is  a  result  of  viscosity. 
Inus  it  seems  inappropriate  to  use  tne  nayleign  criterion  on  tnese  types 
of  flow. 

otability  calculations  for  tne  plane-oscillating  boundary  layer 
(von  Kerszeck  and  Davis  ( 1 97*0 ,  Davis  ( 1 S76 > y ,  or  the  plane  oscillating 
Poiseuille  flow  (Grosch  and  i&lwer.  (i960))  or  tne  round  oscillating  flow 
in  a  pipe  (Yang  and  Yih  (1977))  indicate  that  tne  imposed  pulsations  do 
not  oatce  the  flow  less  stable.  Yang  and  Yih  have  found,  tnat  the  oscil- 
lrtir.7  pipe  flow  is  stable  to  cviniametric  disturDs.nces  at  all  frequen¬ 
cies  and  at  all  heynolds  manners  (basea  in  tnis  case  on  the  amplitude  of 
tne  velocity  oscillations).  no  theoretical  analysis  is  available  for 
tne  stability  of  a  pulsating  pipe  flow,  in  wnich  oscillations  are  super¬ 
imposed  on  a  steady  parabolic  profile.  The  linear  stability  analysis 
applied  separately  to  the  mean  ana  oscillatory  components  of  the  flow  is 
iacapaole  of  predicting  the  outcome  because  the  velocity  profile  enters 
into  urr-Soomerfeld  equation  in  a  nonlinear  way. 

The  first  experimental  investigation  of  velocity  and  pressure  in 
turbulent,  pulsating  pipe  flow  In  water  was  made  by  Schultz-Grunow 
(1940).  The  veloolty  was  measured  by  a  Pitot  tuoe,  and  was  limited  to 
time  averaged  values.  It  was  concluded  that  tne  instantaneous  velocity 
profiles  were  similar  to  the  steady  profiles  in  a  converging  ohannel 
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curing  tae  acceleration  stage  of  a  period,  anc  to  tne  steaoy  profiles, 
existing  in  a  oiversent  cnannel .  it  appeared  tnat  tne  intensity  of  tur¬ 
bulence  increased  curing  tne  aecceleration  stage,  tnis  ooservation  was 
eased  on  large  increase  in  the  scatter  of  tne  Gata  curing  acceleration. 

Cheng  (1971)  used  hot-wire  anemometry  for  velocity  measurments  in 
pulsating  pipe  flow  in  air.  he  nas  found,  tnat  there  is  a  strong  incre¬ 
ase  in  the  Reynolds  stress  and  the  turbulent  kinetic  energy,  resulting 
frox  tne  addition  of  puisations  to  the  steady  flow,  a  strong  nonlinear 
interaction  was  found  to  play  an  important  role  in  the  uistrioution  of 
tne  turoulent  energy,  it  should  be  noted  nowever,  tnat  tne  pipe  usea  in 
Cneng* s  experiment  was  very  short,  about  30  diameters  only,  and  tne  pul¬ 
sations  were  not  simply  naraonic  because  of  tne  awKvard  manner  in  wnicn 
tr.ey  were  introduce^  Consequently  it  is  impossible  to  obtain  reliable 
pne.se  averaged  Information  from  tnese  resuits. 

uerraru  (1971)  investigated  the  turbulent  pulsating  pipe  flow  oy 
flow  visualization  teennique  wnicn  enaded  him  to  ootain  mainly  qualita¬ 
tive  data,  he  found  a  similarity  between  laminar  and  turbulent  oscil¬ 
lating  velocity  profiles.  Uerrard  also  observed  tne  pease  dependence  of 
turbulent  intensity  and  found,  that  all  turbulence  nearly  disappears 
curing  the  acceleration  stage.  The  similarity  oetween  laminar  and  tur¬ 
bulent  oscillating  velooity  profiles  was  noticed  also  by  unamann  (1973), 
wno  used  six  thermistor  prooes  for  instantaneous  velocity  measurments  at 
different  radial  positions. 
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A  niucoer  of  investigators  (oergeev  <  1 9t>t> ; ,  Clarion  and  Peiissier 
(1 9/5),  ner-cii  ana  Thomann  ( 1977 ) ,  among  otners )  tried  to  determine  tne 
transition  Reynolds  numoer  in  oscillating  pipe  flow,  wnenever  the  flow 
oecomes  turoulcnt  during  a  portion  of  the  cycle  only,  a  large  scatter 
in  tne  results  indicates,  that  tne  instability  nas  a  nonlinear  character 
and  depenas  strongly  on  the  experimental  facility,  in  tne  work  of  herk- 
11  and  Thomann,  for  example,  a  pipe  which  was  closed  from  both  sides  wa3 
used;  that  resulted  in  tne  appearance  of  shock  waves  and  the  pressure 
gradient  was  not  spatially  unifona.  Their  results  seem  to  be  relevant 
to  tne  specific  experimental  set-up  only.  All  the  aoove  mentioned  ex¬ 
periments  show  that  the  flow  can  oe  turbulent  during  some  part  of  tne 
cycle,  and  laminar  during  another  part  of  the  cycle. 

tlamen  and  hinton  (1977)  investigated  experimentally  oscillating 
ana  pulsating  pipe  flow  in  water  by  a  hydrogen-ouoole  teennique.  A  good 
agreement  with  theory  was  observed  in  laminar  flow  at  low  Reynolds 
number;  at  higher  Reynolds  numbers,  their  pipe  wnlch  was  only  170  diam¬ 
eters  in  length  was  far  too  short  for  the  flow  to  oecome  fully  devel¬ 
oped.  They  ooserved  that  tne  interaltteney  of  the  pulsating  flow  at 
higher  Re  (Re*2900)  depended  on  both  the  mean  Re  and  amplitude  of  velo¬ 
city  oscillations. 

wizushlna  et  al  (1973*,  1S73b,  1975)  and  haruyama  (1974)  investi¬ 
gated  a  pulsating  turbulent  flow  by  electrochemical  method,  which  en¬ 
ables  meaaurments  of  instanteneous  velooity.  Tne  measured  data  was  re¬ 
corded  on  a  magnetic  tape  for  subsequent  processing  on  digital  computer. 
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iney  found  a  critical  value  of  tne  pulsation  psriou  i  ,  so  that  for 

or 

slow  pulsations  (i.e.  for  1>1‘  )  there  is  no  significant  change  in  tur- 

cr 

oulent  intensity  at  different  phase  angles ,  anc  the  flow  oenaves  like  a 

"steady"  turbulent  flow,  for  nigner  frequencies  (£<T  J,  tne  turbulent 

c* 

characteristics  are  strongly  dependent  on  the  phase  angle,  k  calculated 
ecicy  viscosity  e*x/ (3 u/3r) attained  negative  values  at  certain  pnase  an¬ 
gles,  when  relaainarisation  took  place  during  tne  acceleration.  An  em¬ 
pirical  relation  for  critical  period  of  pulsations  l  was  determined  to 
oc: 


T  C/0a0.1b«e*/3  1.12 
cr 

vr.erc  U  is  tne  tine  mean  bulk  velooity.  it  is  not  clear  from  tnese  re¬ 
sults  wnetner  tus  change  ir.  the  flow  oenaviour  results  from  tne  differ¬ 
ent  frequencies  of  pulsations,  or  there  Bight  be  an  influence  of  the  am¬ 
plitude  of  tne  pulsations  at  hlgner  frequencies. 

i.irase  ( 1 V7v >  used  laser  Doppler  velooioeter  to  measure  pulsating 
turouient  pips  flow  in  water  at  high  values  of  tne  diaensionless  fre¬ 
quency  parameter  (55<  o<137).  Ihe  phase  shift  between  pressure  and  vel¬ 
ocities  at  all  radial  positions  and  at  all  frequencies  was  found  to  be 
SO0.  An  ebay  viscosity  aodel  was  used  to  oaiouista  tne  tiae-depenuent 
velocities  in  this  type  of  flow,  and  a  reasonable  agreement  was  ootaineu 
for  velocity  profiles  at  two  pnase  angles. 

naaaprian  and  snuen-kel  Tu  (ISoO)  ooserved  laminar! cation  of  ini- 
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tialiy  turoulent  flow  of  oil  In  s  pipe  wnen  pulsations  were  Intro-uccv. . 
x.icy  ooser/ec  tr.at  t  **C  max  lx  M wL»i  w  mplituuo  of  tne  velocity  pulsations  near 
tne  surface  of  tne  pipe  is  attained  oefore  tne  maximum  amplitude  occurs 
ir.  t*:o  central  region  of  tne  pipe  wnenever  tne  flow  is  fully  turoulent. 
foor  angular  resolution  of  tne  data  in  their  experiment  dia  not  allow 
tnem  to  get  quantitative  information  aoout  phase  shift  angles.  The  vel¬ 
ocity  profiles,  measured  in  laminar  pulsating  flow,  were  In  good  agree¬ 
ment  with  the  tneoretioal  results  of  Ooniaa.  similarity  was  notes 
oetween  laminar  and  turbulent  flows  at  tne  same  frequency,  and  it  was 
concluded  that  tne  imposed  oscillations  nave  no  effect  on  tne  time  mean 
properties  of  tne  flow. 


i*a“c  It. 
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2. 1  General  Description  of  tne  Experimental  Apparatus 


A  straight  and  smooth  aluminium  pipe  33  mm  in  diaoeter  and  17W1G'3 
am  long  was  used.  The  facility!  snown  schematically  on  fig. 2.1,  was  or¬ 
iginally  used  by  hygnanskl  and  Champagne  (1973),  and  described  in  detail 
in  their  paper.  The  pipe  was  carefully  alligned  to  within  1  mm  over 
it’s  entire  length.  The  contraction  was  made  of  two  subsections,  giving 
an  overall  area  ratio  of  .340:1.  As  a  result  of  the  careful  alignment 
and  tne  smooth  inlet,  laminar  flow  could  be  retained  at  the  Reynolds 
numoers  exceeding  20000  without  the  addition  of  soreens.  The  mean  flow 
was  supplied  by  a  high  pressure  source  (6  atm  compressor;  and  controlled 
by  a  precise  pressure  regulator.  This  arrangement  insured  that  that  tne 
flow  rate  was  independent  of  the  superimposed  pulsations  and  the  flow 
regime ,  wnetner  laminar  or  tubulent,  in  tne  pipe. 

Pressure  oscillations  were  introduced  by  a  valveless  piston  pump, 

connected  to  the  settling  ohamber.  The  piston  diameter  Dp  was  90  mm, 

and  it's  displacement  r  could  be  changed  from  about  r  »5  no  to  r  *75  ma 

P  p  p 

in  17  steps.  The  length  of  the  ecotoh-yoke  was  *p«300  an.  The  ouIk 
rate  of  ohaige  of  the  settling  chamber  for  the  olroular  frequency  w  is 
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it-  c 

.•iver.  oy 

ov/cit =  w  i/'Ar  <»)(ooswt-r  /iA  sini.wtj 
P  P  P 

inc  amplitude  ratio  In  tne  oulit  displacement  octvcen  tne  fundamental 
frequency  w  arc.  It' a  first  narmonic  2to  is  therefore  rp/lA^<o >  at  tne  ni- 
cnest  possible  displacement  amplitude.  The  oscillations  at  tne  namonic 
frequency  could  tnus  oe  neglected.  The  pucap  was  driven  oy  1.5  np  vari¬ 
able  speea  motor,  permitting  a  cnange  in  tne  period  of  pulsations 
octw6en  0.5  see  to  5  sec.  ins  repeatability  oi  tr.e  period  was  Detter 
than  O.J>. 

c.L  i'ne  riston  Pump  as  a  oouroe  of  Pressure  an-  Velocity  Pulsations 

ine  velocity  one  pressure  pulsations  lr.  tne  pipe,  resulting  from 
tne  ao/cment  of  tr.e  piston,  oan  oe  easily  analysed  by  assuming  tnat  tne 
flo*  in  tne  pipe  is  laminar  and  fully  developer.  Tne  amplitude  of  tne 
pulsations  o*  tr.e  ouIk  velocity  is  rhere  n  is  tne  raoius  of 

tne  pipe,  anc.  is  tne  amplitude  of  tne  oscillations  in  tne  flow  rate. 
Uj  is  proportional  to  dPj/dx,  wnere  p1  is  the  amplitude  of  pressure  pul¬ 
sations  at  tne  entranoe  of  the  pipe.  for  fully  developed  flow 
Sp^/axsp^/L,  vnere  L  is  the  lengtn  of  tne  pipe,  for  a  given  amplitude 
of  pressure  oscillations  the  amplitude  of  tne  ouIk  velocity  Uj  depends 
strongly  on  frequenoy.  If  T7  la  tne  steady  part  oi  tne  volumetric  vcio- 
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city  {Csfy  ffr.**,  v-r.'rc  C  is  tne  mean  iiov  rate;  a  no  p  is  tne  steaay  pres¬ 
sure,  tncr.  tne  relationship  cctvc.n  tne  ensamcie  wverageu  time  aepenaent 
parts  <c>  E.nu  <p>  can  oe  expresses  oy  tne  equation 

<U>/U=<d>/».;  =  <J  i<P>/pjCOS  4>  2.1 

Q 

where  the  amplitude  coefficient  a  s(.*t./Q)/(p ,/p;  ano  pnasc  lag  angle  A 

q  1  i  s 

were  calculates  for  tne  fully  dcvelopeu  laminar  flow  oy  Ucniaa  (Isco) 
are  reproducea  in  rig.  2. 2.  tne  aoscissr.  in  tnis  Figure  is 
lop  o=iogK*Wv.  it  can  he  seen  from  Fig .2. 2,  tnat  at  very  low  frequen¬ 
cies  {  a<1)o  *1  ano  tne  pnasc  lag  6  -K).  ’inu3,  wnen  tne  frequency  of 

<i  <5 

tne  imposed  pulsations  is  low,  the  flow  at  any  instant  Denaves  litce  a 
Foiseullle  flow  at  the  appropriate  instanteneous  pressure  gradient. 
Inertia  effects  oecoae  noticeable  with  increasing  frequency  wnen  the 
flow  cannot  follow  after  the  rapia  changes  in  pressure  anymore.  iJi  in¬ 
crease  in  a  causes  a  corresponding  increase  in  tne  prase  angle  oetween 
tne  volumetric  velocity  and  tne  pressure  <J>  ,  and  dampens  tne  amplitude 
of  tns  veiocity  pulsations,  as  well,  at  very  nigr.  frequencies  (a  »1) 
the  inertial  term  becomes  dominant  In  comparison  witn  tne  viscous  term, 
ana  tne  resulting  acceleration  (dQ/dtsiu)  Gj)  is  proportional  to  tne 
driving  force  ( instanteneous  pressure).  The  quantity  iGj  is  tnus  in 
pnase  witn  tne  pressure,  wnile  the  pulsations  in  flow  rate  nave  *0° 
phase  lag  benina  p  ^ . 

consider  tne  control  volume  in  wnicn  tns  movement  of  the  piston 
causes  a  corresponding  change  in  tne  volume  of  the  settling  cnamoer. 
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if  <U>  is  tine  time  dependent  component  of  the  volume  flow,  tne  continui¬ 
ty  equation  for  tne  control  volume  may  oe  expressed  oy: 

dn/dtsd/dt(  p V)=-<Q  p >  2.2 

where  V  is  tne  volume,  pis  tne  density,  and  <Q>  is  assumed  to  oe  posi¬ 
tive  when  the  fluid  leaves  the  settling  chamber  tnrougn  tns  pipe. 

for  small  amplitudes  of  volume  oscillations,  (relative  to  the 
total  volume  of  the  settling  chamber  VQ)  ana  pressure  oscillations,  Pj, 
(relative  to  tne  ataospneric  pressure  po>  equation  (2.2)  can  oe  linear¬ 
ized  to  give 

p  dV/dt+V  d  (/dt=-p  <Q>  2.3 

o  o  o 


equation  (2.3)  is  a  linear  homogeneous  equation,  naving  solutions 
of  the  form 


V«Vo+Vfexp(i(  wt*  <>v)) 

Z.H 

P»P0^P*P?«P<i“  t) 

2.0 

,  «5>*Q,exp(i(  ut-  6  ) 

1  q 

2.0 

i 


‘ .at;  i. 
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where  <Pv  and  represent  the  phases  of  the  volume  and  flow  pulsa¬ 
tions,  respectively,  relative  to  the  phase  of  the  pressure  osoillation 
which  was  arbitrarily  set  to  be  zero.  Assuming  that  air  obeys  the  ideal 
gas  law 

prp^pX  2.7 

where  Ik  is  gas  constant,  and  that  the  process  is  isothermal,  and  substi¬ 
tuting  equations  (2.4)  to  (2.7)  into  (2.3)  one  oo tains 

PQiu)  Vfexp(i(u  t+$  v)+Vo1oj  PjexpUui  t)+pQQ1exp( i(  ojt-  <J>(J))=0  2.6 

Eq .(2.6)  contains  two  unknown  amplitudes,  p1  and  ,  but  the  as¬ 
sumption  of  fully  developed  laminar  flow  provides  an  additional  equation 
relating  the  two  quantities.  Substituting  (2.1)  into  (2.6)  and  cancel¬ 
ling  exp(iut)  yields 

iw  poVjexp(i  $y)«-p1(iu>  VQ-  1/p  SpQ  oqexp(-i  $q))»0  2.9 

Equation  (2.9)  can  be  separated  into  real  and  imaginary  parts  and  solved 

H 

for  $v  and  p^  In  laminar  pipe  flow  0/p»  irD  /126  v>L  is  a  known  quan¬ 
tity  fixed  by  the  flow  geometry  and  the  viscosity  of  the  fluid. 
Substituting  it  into  (2.9)  and  taking  real  part  of  the  equation  gives 

-  wVjSin  ♦y+(Q/p)P1  oqcos  $q«0 


2.10 
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ii”  i  in.-ry  part  oi  (2.*;  provides  tne  secern  equation 


cos  <J>y 


Qo  sin  A 

q  _  q  j 

pw 
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One  -ay  now  cxspress  4> y  ano  in  terms  of  the  piston  displacement 
ano  frequency  oi  pulsations  ui 


tg  K 


P  ^  o  cos  6 
o  -  q  Yq 

_ p_ 


pi  = 


(u)V  +  P  ^  cr  sin  4>  ) 
o  o  =  q  q 

uV1  sin  $v 


P  q 


a  cos  4> 


2.12 
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ine  amplitude  of  tne  oscillations  i.i  the  flow  rate  Q1  can  oe  expressed 
in  terms  of  p1  an^  hq.(2.1) 


Ct1  =  <P1/p)CT(.C 


■.14 


for  a  given  amplitude  of  the  volumetric  change  of  tne  system  V^, 
amplitudes  of  tne  oscillations  in  tne  flow  rate  tj  ana  in  pressure  p1 
were  calculated  from  equations  (2.12)  to  (2.14)  as  functions  of  the  per¬ 
iod  (fig.2.3a).  at  iow  frequncies  (large  values  of  the  period  'i,a  <1) 
cr  ~1,  and  the  pressure  amplitude  in  the  settling  cnamoer  is  adequate  to 
"pusn"  tnrougn  the  pipe  nearly  all  tne  mass  flow  supplied  by  tne  piston 
pump  to  tne  settling  chamber.  At  these  frequencies  tne  amplitude  in 
flow  rate  is  tnus  proportional  to  tne  rate  of  change  in  tne  volume  of 
tne  settling  chamber  V^.  for  <*<1  therefore  tne  relation  1/T  is 
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obtained. 

At  higher  frequencies  (a  >1)  decreases  witn  increasing  frequency 
(see  Pig. 2. 2),  and  the  amplitude  of  the  pressure  oscillations  in  the 
settling  chamber,  resulting  from  the  oscillations  in  the  volume  of  the 
system  grow  more  steeply  than  the  flow  rate.  In  order  to  understand  the 
physical  meaning  of  the  mutual  interaction  between  the  oscillations  in 
the  pressure  and  in  the  flow  rate,  equation  (2.6)  is  rewritten  in  vector 
form,  using  the  notation  ^ ex p C  14> y) ;  and  Q=pQQ1exp(-i<|iq) : 

iV+iP+QrO  2.6a 


The  calculations  indicate,  tnat  the  oscillations  in  pressure  and  volume 
are  nearly  160°  out  of  phase  (<{iv%1o0o),  i.e.  pressure  increases  when 
the  piston  pushes  air  into  the  settling  chamber,  and  vice  versa.  At 
frequencies  corresponding  to  as  5  the  resulting  oscillations  in  the  flow 
lag  behind  the  pressure  by  more  than  70°  (see  Fig. 2. 2).  Tawing  into  ac¬ 
count  these  phase  relations,  and  the  fact,  that  multiplication  by  i  adds 
90°  to  the  phase  angle,  two  sketches  can  be  drawn,  the  first  showing  the 
relative  direction  of  each  vector,  and  the  second  representing  the  mass 
balance  in  accordance  with  Eq.(2.tta): 


V 
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it  follows  from  tne  diagram,  tact  tne  pressure  vector  iP  is  in- 

-►  Q 

clinea  to  tne  flow  rate  vector  s,  oy  aoout  lou  .  mien  tae  pressure  in 
tne  settling  cnamoer  is  maximum,  tne  time  cepcnccnt  portion  of  tne  flow 
rate  is  directeo  towards  tne  settling  cnamoer,  causing  a  furtner  incre¬ 
ase  in  pressure  amplitude,  lne  pnase  relations  oetween  pressure  ana 
flow  rate  tnus  lead  to  a  wind  of  a  positive  feeaoscK,  resulting  in 
curves  Q^sQ^(  a)  and  p^sp  (a  )  (fig. 2oa>,  which  resemDle  a  resonance 
pnenomenon.  The  amplification  of  pressure  oscillations  nas  however 
nothing  in  common  with  resonance,  since  in  tnis  case  tne  system  has  no 
eigen  frequency. 

at  even  higher  frequencies  (a  >o)  is  small,  ana  tne  pressure  am¬ 
plitude  is  insufficient  to  prcvice  strong  oscillations  in  the  flow  rate, 
tnus  Cij  decreases,  and  tne  amplification  influence  of  tne  flow  rate  on 
tne  pressure  fluctuations  weakens;  so  the  pressure  amplitude  decreases 
togetner  witn  the  amplitude  of  flow  rate,  furtner  increase  in  frequency 
(a  >10)  results  in  a  0,  and  vanishing  oscillations  in  flow  rate.  In 
tnis  case  tne  entrance  from  tne  settling  chamoer  to  pipe  may  be  regarded 
as  closed  for  the  time  dependent  flow,  and  tne  amplitude  of  the  pressure 

oscillations  tends  to  a  constant  value  p  V,/V  wnicn  correspond  to  the 

ol  o 

ideal  gas  law  pVscorst. 

The  results  of  Fig.  2.3a  were  calculated  for  tne  geometrical  data 
appropriate  to  tne  existing  experimental  facility.  Tne  total  volume  of 
tne  settling  chamber  was  approximately  equal  to  1  m3,  and  the  maximal 
changes  in  volume  were  100  cm13.  The  preceeding  calculations  show,  that 
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the  maximum  responce  of  tne  system  to  the  volume  oscillations  occur  at 
T=5.2o  sec  (a=-_>.7)  (Fig.2.3a) .  in  reality,  the  maximum  amplitude  of 
pressure  pulsations  occured  at  1=2.4  sec,  corresponding  to  a  sb.7. 
(Fig. 2. 3°).  The  discrepancy  between  theory  and  experiment  is  attributed 
to  the  Influence  of  the  entrance  region  which  was  not  considered  in  the 
calculations,  and  some  difference  in  the  geometrical  parameters. 
Fig. 2. 3  shows  that  the  calculated  maximum  amplitude  of  the  pressure  os¬ 
cillations  is  aoout  three  times  higher,  than  it  woulu  nave  been  if  the 
settling  cnamber  was  closed  (i.e.  A  p(  a=5.7)/  Ap(a-*®°  )  :j)  The  maximum 
amplification  in  pressure  oscillations  realized  in  laminar  flow  was  ap¬ 
proximately  2.5,  and  the  absolute  value  of  Ap  at  a=o.7  agrees  fairly 
well  with  tne  calculations  (see  Fig.2.3b). 

borne  amplification  exists  in  turbulent  pulsating  flow  also,  as  it 
can  be  seen  from  the  Fig .2, 3b,  but  the  maximum  amplitude  is  reduced. 
The  mechanism  of  the  amplification  in  turbulent  flow  is  qualitatively 
identical  to  the  mechanism  in  laminar  flow.  It  will  be  snown  in  Section 
5.2,  the  oscillating  parts  of  the  bulk  velocity  <u>  in  ooth  laminar  and 
turoulent  flow  behave  alike,  but  the  phase  lag  angle  is  notably  less 
in  the  turbulent  flow,  than  in  the  corresponding  laminar  case.  The  com- 
ponent  of  the  flow  rate  vector  Q,  wnioh  is  coll inear  with  the  pressure 
vector  IP  and  thus  causes  the  amplification  (see  the  sketches),  has  the 

length  |Q|sln6  ;  and  thus  a  lower  value  of  <t>  in  tne  turbulent  case 

<!  Q 

causes  a  reduction  in  pressure  amplification.  The  maximum  of  the  curve 
p^sp^a)  is  shifted  in  «  relative  to  its  position  in  tne  laminar  flow 
obviously  due  to  different  amplitude  a  «o  (a)  and  phase  d>  *<fc  (a)  rela- 

Q  Q  Q  Q 
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tions  in  this  case. 


2.5  The  Influence  of  tne  Lengtn  of  the  Pipe  on  the  Pulsating  Flow 


It  is  usually  assumed  tnat  in  a  fully  developed  region  of  a  pulsat¬ 
ing  pipe  flow  the  pressure  gradient  depends  on  time  only,  and  conse¬ 
quently  the  radial  component  of  the  pulsating  velocity  has  to  vanisn. 
Only  recently  Richardson  (19o0)  raised  this  question  in  his  discussion 
of  the  paper  of  Kirmse  (1979)-  Ricnardson  argued  that  applying  a  pul¬ 
sating  pressure  at  one  end  of  the  pipe  does  not  necessarily  result  in  a 
linear  pressure  distribution  along  the  pipe;  acoustic  waves  may  make 
pressure  gradient  dependent  on  tne  axial  coordinate  and  resulting  in  the 
generation  of  a  radial  component  of  velocity.  It  was  noted  in  Section 
1.4,  that  the  fully  developed,  laminar  pulsating  flow  is  independent  of 
the  axial  coordinate  because  the  inertial  term  in  the  Navier-Stokes 
equations  (1.4)  (uV  )u  dissapears,  making  the  equation  linear  and  ena¬ 
bling  a  separation  of  the  steady  from  the  oscillating  flow  components. 
Tne  pressure  gradient  is  independent  of  the  axial  coordinate  whenever 
any  change  in  pressure  at  the  inlet  of  the  pipe  is  felt  "Instantly”  in 
tne  entire  pipe.  Since  a  weak  pressure  pulse  travels  at  the  speed  of 
sound,  a,  the  assumption  of  spatially  constant  pressure  gradient  may 
only  oe  valid  for  a  limited  length  of  the  pipe  L. 


In  order  to  obtain  an  estimate  of  the  highest  permissible  pulsation 
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frequency  and  tne  longest  pipe  length  L  for  which  the  pressure  drop 
3p/  3x  may  be  assumed  constant  tne  flow  is  considered  to  oe  incompres¬ 
sible  (i.e.0«a).  The  pressure  is  prescribed  at  ootn  ends  of  the  pipe: 
psO  at  tne  open  ena  of  the  pipe  at  x=L;  and  psp+p^osO*)  t)  at  the  inlet 
to  the  pipe  at  x=0.  When  U«a  and  can  be  neglected,  one  obtains  the 
wave  equation  for  tne  spatial  and  temporal  pressure  distribution 

o‘‘p/dt2=atdtp/dxt  2.14 

The  solution  of  (2.14)  for  tne  i^iven  boundary  conditions  is: 
p  sm  ~  (L-x) 

p(x,t)  =  -  (L-xJ+Pj  - - - sintut  2.15 

L  sin  ■-  L 

a 

Pressure  gradient  tnerefore  becomes: 


iE 

dx 


(x.t) 


0,  cos  a  «-»> 

a  sin  L 
a 


sinwt 
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Tne  pressure  gradient  may  be  regarded  as  spatially  constant  provided 
x  o)/a«1  for  all  x  concerned,  or  tnat  tne  lentil  of  tne  pipe  L  is  mucn 
snorter  than  the  acoustic  wave  lentil  L<<2  Tia/ui  . 


A  more  accurate  estimate  of  the  limitation  on  tne  frequency  of  os¬ 
cillations  of  a  slightly  compressible  fluid  in  a  pipe  of  a  given  lengtn 
is  obtained  oy  considering  the  continuity  equation  for  the  pulsating 
part  of  tne  bulk  velocity  <u> 
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ft  (PQ^<P>)  ♦  (PQ+<P>)  = 
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vnere  <  t  y  is  tne  oscillating  part  of  tne  wall  shear  stress 
<  tv>*-pl  3<u>/  3r)w.  ey  assuning  tnat  tne  fiuio  is  only  siigntiy  com¬ 
pressible,  (  i.e.  <p>/  pQ«1>,  equations  (2.17/  anc  ^2.1b)  can  be  line¬ 
ar  izeu  to  give 
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me  pressure  tern  can  be  eliminated  froc  (t.20>  using  the  adiaoatlc 
sound  velocity 


me  wall  snear  stress  <  t^>  may  be  eliminated  oy  following  the 


equation 
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and  equation  (2.1),  where  and ~p  are  the  steady  wall  shear  stress  and 
pressure,  respectively,  bhear  stress  and  pressure  pulsations  amplitude 
ratio  and  relative  phase  shift  ^  were  calculated  using  the  anallti- 
cai  solution  of  Uchida  for  the  laminar  pulsating  flow.  From  tne  equa¬ 
tions  (2.1)  and  (2.22) 


y  °t  ui  ,i(w 

\  °q  0 

The  steady  bulk  velocity  U  and  tne  steady  wail  shear  stress  t  are  re¬ 
lated  by  the  equation 
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Re 


He  therefore  obtain 


TW*°  PQ  V/R  /  0q)eXP^1^x  + 


By  assuming  harmonic  oscillations,  the  time  derivatives  can  be  re¬ 
placed  by  the  operator  iw  ,  so  that  3<p>/3tsiwp  and  3<U>/3t*iw  u^. 
Equations  (2.19)  and  (2.20)  were  therefore  rewritten  in  the  form: 


2.24 


and 


dpY  j  8p  v  o  !(♦♦♦) 

dr--7(1*,’.*Trre  >ui 

a  a 


c. t.t).' 


•‘‘V-'-iM 


-  ci/-JC  =tnpi 


*1*  1U '  pc 


iuip  a  e(4  +<fc  ) 

Vy(1+!rre  q> 

a  a 


tc  .  C.  4.  . 


-utrooucin*  Ys^i,  trie  soiutior.  ol  {>.. 


•  t  err  oc  prrsci.tCi.  ii.  t:.c  ic-rv 


b^cA^cxpi-  Y:  ;+i-cxp(  Yx) 


. sin?.  (L.tb)  er.u  U.c.o;  tus  censity  distribution  oe coses 


u?i-  Yx)-i.,/i  aqi  Y") 

110  c  0 


v-nerr.  i-  =  t./i,. 
o  1  1 


•  ■*  **<« 


Page  30 


'io  determine  the  integration  constants  and  the  boundary  con¬ 
ditions  at  both  ends  of  the  pipe  are  used.  At  the  outiet  the  amplitude 
of  tne  pressure  pulsations  vanishes  and  so  does  thus  the  amplitude  of 
density  pulsations  P  (xaL)sO.  for  the  given  amplitude  of  pulsations  in 

tne  flow  rate,  C^,  the  boundary  condition  on  the  amplitude  of  the  bulk 

2 

velocity  pulsations  at  the  inlet  is:  U^(xsO)ag  /  hr  .  Taking  into  ac¬ 
count  the  boundary  conditions  the  final  form  of  the  solution  is: 


U  -  _1  cosh  y(L- 
1  "  7TR2  cosh  r 

p  =  sinh  Y(L-: 

1  ttr2  cosh  n 


2.30 


2.31 


As  a  quantitative  measure  of  the  flow  dependence  on  frequency  a 
phase  shift  between  tne  bulk  velocity  pulsations  at  both  ends  of  the 
pipe  was  chosen.  It  is  obvious  that  for  the  flow  to  be  Independent  of  x 
this  phase  shift  has  to  vanish.  The  complex  parameter  which  governs  the 
phase  of  U  at  different  locations  along  the  pipe  is  Y ,  which  in  itself 
depends  both  on  the  length  of  acoustic  wave  A=2tra to  and  on  the  fre¬ 
quency  parameter  a .  The  parameter  is  dependent  on  the  radius  of  the 
pipe  and  on  the  viscosity  of  the  fluid.  The  major  factor  governing  the 
phase  shift  was  found  to  be  the  ratio  A  /L  (i.e.  the  aooustio  wave 
length  to  the  lengtn  of  the  pipe),  fig. 2. 4  gives  numerioal  results  of 

solution  (2.30)  for  the  parameters  used  in  this  experiment  (IU1.65 

2 

cm, 1*1700  cm)  and  v*0.16  cm  /sec.  It  can  be  seen  from  fig. 2. 4  that  for 
L/A<0.08  the  resulting  phase  shift  is  less  than  1*  .  Accepting  this 
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value  as  the  upper  limit  for  whicn  one  may  assume  3p/  3x=const.  gives 
the  maximum  permittee)  frequency  parameter  a=15,  and  pulsation  period 
longer  than  0.5  see.  It  is  worth  noting  tnat  for  L/X  >0.2  the  flow  can 
no  longer  be  regarded  as  one-dimensional,  even  by  a  crude  approximation. 


2.4  Calibration  of  Velocity  and  Pressure  Sensors 


Velocity  measurments  were  made  with  a  rake  of  9  hot  wires,  distri¬ 
buted  evenly  in  the  radial  direction  at  distances  equivalent  to 
&r/R=0.12  between  the  neighbouring  wires;  so  that  when  the  first  wire 
was  on  the  center  line  of  the  pipe,  the  y-th  wire  was  located  at  a  dis¬ 
tance  0.5  mm  from  the  wall  (i.e.  at  r/R=0.97).  All  velocity  measur¬ 
ments  were  taken  at  the  exit  plane  of  the  pipe.  A  10  channel  constant 
temperature  hot-wire  anemometer  and  10  channel  amplifier,  both  built  by 
the  electronic  shop  of  the  School  of  Engineering  of  Tel-Aviv  University 
were  used  in  the  experiment.  The  outputs  of  the  amplifier  were  connect¬ 
ed  via  an  analog  to  digital  converter  to  a  DEC  POP  11/60  minicomputer. 

The  calibration  of  the  hot  wires  is  done  in  a  wind  tunnel,  which 
provided  a  stable  velocity  stream  between  30  cm/sec  and  15  m/sec. 
Although  flow  reversal  was  avoided  in  the  experiment,  very  low  veloci¬ 
ties  occured  in  the  pipe  as  a  result  of  the  superimposed  pulsations.  It 
was  thus  necessary  to  calibrate  the  wires  at  the  lowest  velocities  anti¬ 
cipated  in  the  experiment. 
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It  is  well  renown,  nowever,  tnat  standard  Pitot  tuoe  is  not  a  proper 
instrument  for  measuring  air  velocities  oeiow  I.p  m/s.  An  alternative 
way  for  accurate  measurment  of  low  velocities  may  be  cased  on  the  fre¬ 
quency  of  vortices  sned  behind  a  circular  cylinder.  Insert  in  Fig.  2.6 
(Kovazsnay  (19**y))  3hows  the  dependence  of  the  frequency  of  shedding, 
i.e.  the  Strouhal  number  St=fD/U,  on  the  Reynolds  number  of  the  cyl¬ 
inder  suggesting  that  St  is  independent  of  He  for  He>300.  The  later 
measurments  by  Roshxo  (1956,  IS0 1 )  confirmed  the  results  snown  in 
Fig. 2. 5  for  he>300.  At  Re<300  St  decreases  with  decreasing  He;  and  the 
scatter  in  the  data  collected  by  various  investigators  is  very  large  in 
this  area.  Different  empirical  formulae  showing  the  dependence  St=f(He) 
(see,  for  example,  Goldstein  (19b5),  Berger  and  wille  (1972))  are  not 
very  accurate  and  may  cause  significant  errors  at  low  Reynolds  numbers. 

In  this  experiment  a  cylinder  1.269  cm  in  diameter  was  used  for  tne 
velocities  under  consideration  (U>35  cm/sec);  He  thus  was  larger  than 
300  and  Strouhal  number  therefore  could  be  regarded  as  a  constant.  The 
dependendence  of  St  on  He  was  measured  for  Re>1700,  where  a  Pitot  tube 
could  provide  accurate  measurments  of  velocity,  and  it  was  assumed  that 
the  mean  value  of  St  could  be  extrapolated  to  lower  Reynolds  numbers. 

A  hot  wire  was  placed  several  diameters  downstream  of  a  cylinder  at 
a  radial  location  at  which  the  vortex  shedding  frequency  is  detected 
most  clearly  on  an  oscilloscope,  and  the  output  of  the  anemometer  was 
sampled  digitally  by  the  computer.  The  sampling  frequency  was  fixed  by 
the  signal  frequency,  which  was  roughly  estimated  from  monitoring  tne 
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oscilloscope.  The  minimum  sampling  frequency  chosen  was  approximately 
twice  the  expected  signal  frequency.  A  Puffer  of  2Q4o  sampled  data  po¬ 
ints  was  Fourier  transformed,  and  the  power  spectrum  of  the  signal  was 
calculated.  The  procedure  was  repeated  several  times  (usually  4  to  10, 
depending  on  sampling  frequency),  the  average  specrum  was  calculated, 
and  the  frequency  of  the  most  energetic  component  was  assumed  to  be  the 
signal  frequency.  The  repeatability  of  the  frequency  measurment  was 
better  than  0.2fc  which  seems  equivalent  to  the  resolution  of  the  metnod. 

On  the  Fig. 2. 5  the  measured  dependence  of  St  on  he  is  shown; 
Fig. 2. 6  gives  the  calibration  curve  for  the  dependence  of  flow  velocity 
on  the  vortex  shedding  frequency. 

The  detailed  description  of  the  hot  wire  calibration  procedure  may 
be  found  in  the  thesis  of  Oster  0980)  and  in  wygnanski  and  Oster 
(ISbl).  Seven  calibration  velocities,  determined  by  the  vortex  shedding 
method,  were  used.  A  4-th  order  polinooial,  giving  the  dependence  of 
the  flow  velocity  u  on  the  output  voltage  E  was  found  by  a  least  square 
method  from  the  7  measured  points: 

4  3 

usa^E  *• • • •♦ag 

The  calibration  curve  was  further  cheoked  for  several  additional  veloci¬ 
ties. 


In  order  to  calibrate  an  x-wlre  the  output  voltage  was  sampled  at 
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11  different  angles  (in  the  range  of  27°>B  >-27°)  at  eacn  one  of  7  velo¬ 
cities  used.  The  output  of  each  wire  depends  on  botn  aosoiute  value  of 
tne  velocity  vector  U  and  its  angle  :  E ^  =t ^ (IuX 0  )  and  E2=E2(|ull3  ). 
It  was  assumed  that  from  the  measured  values  of  and  U  and  B  can 
be  uniquely  determined.  From  the  measured  77  calibration  points  the 
4-th  order  polinomials,  giving  the  dependence  of  the  aosoiute  value  of 
the  velocity  vector  juf  and  its  angle  p  on  the  output  voltages  of  both 
wires  Ej  and  E^  were  found 

(u|»a  t  E ,  4+a2E24+ . . . +a ,  3E , +a  1 4£2+a ,  5 

^*b1E11,+‘52B2i*+‘  *  '+b131S1+bl4E2','b15 

With  fu|  and  3  Knowr^  the  velocity  components  were  calculated  from  the 
equations  usfu|co3$  and  v=ju|sin$  . 

A  Valldyne  model  DP215-30  pressure  transducer  with  OP 15  Sine  wave 
Carrier  Demodulator  was  used  to  measure  the  pressure.  The  pressure 
transducer  was  connected  by  a  5  cm  long  tygon  tube  to  tne  second  section 
of  bell-shaped  nozzle.  For  the  calibration  of  the  Demodulator  output  a 
Fuess  micromanoaeter  was  connected  in  parallel  with  the  transducer.  The 
pressure  transducer  provides  linear  responae  in  the  whole  range  of  pres¬ 
sures.  The  dynamlo  responoe  of  the  pressure  transducer  (1000  hz)  was 
more  than  adequate  for  the  present  purpose. 
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2.5  Data  Aquiaition 


A  D.C.  voltage  corresponding  to  the  output  of  the  anemometer  in 
absence  of  flow  was  subtracted  from  the  signal  before  amplification  in 
order  to  take  advantage  of  the  full  range  of  the  A/D  converter  which  ac¬ 
cepts  signals  between  ^5  volt.  The  converter  having  a  12-bit  presision 
provides  a  resolution  of  about  2.5  mv.  The  amplified  anemometer  signals 
together  with  the  output  from  the  pressure  transduser  were  sampled  at  a 
predetermined  frequency,  the  sampled  data  was  converted  into  16-bit 
words  and  arranged  in  buffers. 

At  the  initiation  of  a  measurment  the  period  of  pulsations  was  de¬ 
termined  by  the  computer.  Two  different  methods  to  measure  the  period 
were  used  in  this  work,  the  first  based  on  the  50  Hz  clock  of  the  com¬ 
puter,  and  the  second  on  1  MHz  clock.  An  optical  switch  (Monsanta  MCA8) 
was  used  to  obtain  trigger  signal.  A  cylinder  1.5  mm  in  diameter,  con¬ 
nected  to  the  driving  motor  of  the  piston  pump,  passed  at  each  revolu¬ 
tion  trough  a  narrow  gap  of  the  optical  switch,  causing  change  in  the 
output  current,  which  operated  TTL  Schmidt  trigger.  In  the  first  method 
the  time  ellapsed  between  two  neighbouring  trigger  signals,  supplied  to 
the  interrupt  input  of  the  computer,  was  measured  by  50  Hz  clock.  The 
resolution  of  a  single  measurment  is  20  msec,  the  final  result  was  obta¬ 
ined  by  averaging  10  measurments.  In  the  second  method  the  trigger  sig¬ 
nal  served  as  input  to  the  A/D  converter  and  was  sampled  at  a  predeter¬ 
mined  sampling  rate,  controled  by  the  1  MHz  clock.  The  period  was  der- 
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i .'•rv.  irora  tne  numoer  or  sampled  points  between  two  consequent  trigger 
.  u».ci  t..e  time  interval  between  samples .  The  resolution  of  this 
...  t  ,'-ri  .sp.-nhs  on  sampling  frequency  ar.e  was  in  principle  orders  of  mag¬ 
nitude  ..lore  accurate  tnan  in  tne  first  method .  In  tnls  case  too,  tne 
final  perioe  was  determined  by  averaging  over  10  cycles. 

with  the  period  of  pulsations  Known,  tne  sampiing  frequency  was 
fixed  so  tnat  1024,  2040  or  40$6  points  were  sampled  per  channel  per 
period  in  order  to  facilitate  tne  processing  of  aata  using  a  fast  four- 
ier  iransfora  ( ffT) .  an  estimate  of  the  repitability  of  the  period  from 
one  event  to  another  was  made  counting  tne  number  of  sampled  points  over 
many  cycles,  and  it  was  found  to  oe  better  tnan  0.3>>- 

iwo  different  sampling  methods  were  used.  In  the  first  methoa  a 
trigger  signal,  supplied  to  the  interrupt  of  the  computer,  initiated  tne 
sampling  wnicn  lasted  precisely  for  1  period,  from  tne  sampled  data  po¬ 
ints  appropriate  velocities  and  presssure  were  calculated  and  tne  re¬ 
sults  were  recorded  on  a  magnetic  tape.  After  tne  data  recording  was 
completed  the  computer  was  ready  to  accept  new  information.  AH  x-wire, 
-no  some  of  normal -wire  raKe  aata  was  acquired  using  this  method,  wnicn 
iv:  s  tan  advantage  that  it  does  not  constrict  tne  sampiing  frequency, 
in  re  are,  nowever,  two  drawbaexs.  firstly,  tne  overall  acquisition 
time  Is  much  longer  than  the  time  of  measurment,  because  the  computer 
requires  additional  time  to  convert  tne  signals  to  velocities  and  pres- 
m>rc.:  and  to  record  the  processed  data  on  tape,  lhe  overall  duration  of 
ttjm  u:..asuraent  increasd  thus  oy  a  factor  of  2  in  tne  case  of  tne  normal 
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wires  and  by  a  factor  of  4  for  the  single  x-array.  The  second  drawback 
stems  from  tne  fact  that  the  duration  of  the  measurments  corresponds  ex¬ 
actly  to  one  period  of  pulsations,  thus  no  spectral  information  can  be 
oDtained  for  frequencies  wnich  are  lower  than  the  frequencies  of  pulsa¬ 
tions. 

A  different  sampling  method  was  used  therefore  in  conjunction  with 
a  rake  of  normal  wires.  The  data  points  were  sampled  continuously  dur¬ 
ing  more  than  6  periods  of  pulsations.  The  memory  was  divided  into  two 
buffers,  and  while  one  buffer  accepted  the  information  sampled,  the  con¬ 
tents  of  the  other  buffer  was  recorded  on  a  magnetic  tape,  thus  provid¬ 
ing  the  possibility  of  essentially  unlimited  in  length  continuous  sam¬ 
pling.  The  output  of  the  optical  switch  was  connected  to  an  additional 
input  channel  and  provided  phase  information.  A  total  of  11  data  chan¬ 
nels  were  thus  sampled:  9  channels  contained  velocity  information,  1 
pressure  and  1  phase  information.  This  method  does  not  have  the  two 
drawbacks  of  the  mentioned  before,  but  it  is  limited  to  a  sampling  fre¬ 
quency  which  could  not  exoeed  1600  Hz.  This  frequency  however  was  quite 
adequate  in  the  range  of  He  considered.  It  should  be  noted  that  in  this 
method  the  "raw"  data  was  recorded  on  the  tape,  and  the  calculations  of 
velocity  and  pressure  were  made  at  a  later  stage. 

In  the  latter  type  of  sampling  single  continuous  record  consists  of 
more  than  6  periods.  In  laminar  flow  7  recordes  were  usually  acquired, 
providing  56  measured  periods  of  pulsations  containing  1024  sampled  po¬ 
ints  in  each  period.  In  turbulent  flow,  the  number  of  records  was  usu- 
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ally  30,  giving  240  periods  witn  tne  number  of  points  in  each  period 
ranging  from  1024  to  40yo,  depending  on  the  duration  of  tne  period,  in 
laminar  flow,  tnorefore,  57344  points  were  sampled  for  each  data  channel 
per  measurment,  while  in  the  turoulent  case  tne  number  of  sampled  points 
per  channel  varied  from  2.5*10^  to  10^. 

2.6  Preliminary  Measurments 

In  most  cases  velocity  aeasunnents  were  taken  in  laminar  and  turou¬ 
lent  flow  regimes  consecutively  while  keeping  all  the  flow  parameters, 
(e.g.  mean  heynolds  number,  period  of  pulsations  and  piston  stroke) 
constant.  before  recording  the  data  on  magnetic  tapes,  preliminary 
measurments  were  made,  in  which  the  data  was  sampled  at  tne  rate  of  1024 
points  per  channel  per  period  for  several  periods.  The  time  mean  velo¬ 
city  was  calculated  for  each  channel  by  averaging,  and  mean  Reynolds 
number  and  flow  rate  were  calculated  by  numerical  integration  using  a 
Simpson  formula.  The  exact  distances  between  tne  wires  were  measured 
with  a  microscope  and  were  used  in  the  integration  procedure.  In  lami¬ 
nar  pipe  flow  tne  center  line  velocity  is  double  the  mean  bulk  velocity. 
This  fact,  coupled  with  independence  of  the  mean  flow  rate  from  the 
state  of  the  flow  (whether  laminar  or  turbulent)  was  used  to  allign  the 
ratce  at  the  exit  of  the  pipe.  The  preliminary  measurments  also  served 
as  a  simple  cneck  of  the  drift  in  the  output  of  the  anemometers. 

It  was  observed  that  the  hot-wire  which  was  placed  at  a  distance  of 
0.5  mm  from  the  wall  overestimated  velocity  values  beoause  of  heat 
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transfer  to  the  metai  wall  whicn  was  absent  In  tne  calibration  condi¬ 
tions  in  the  wind  tunnel.  At  tne  Reynolds  numbers  usea,  tne  velocities 
near  the  wall  were  usually  less  than  0.7  m/sec.  Tnus,  tne  heat  transfer 
wa3  greatly  influenced  by  free  convection,  wnien  is  strongly  dependent 
on  the  flow  geometry.  For  this  reason  the  mean  output  of  tne  last  wire 
was  ignored  in  the  integration  procedure  and  served  only  for  obtaining 
information  on  the  oscillating  and  turbulent  components  of  velocity. 
The  non-slip  condition  at  tne  wall  was  taken  as  an  additional  point  for 
the  purpose  of  integration.  The  hot  wire  placed  on  the  center  line  of 
the  pipe  did  not  contribute  to  the  flow  rate  oecause  of  its  vanisning 
radial  distance. 

For  all  radial -raxe  data,  only  seven  measured  velocities  were  used 
to  determine  the  mean  flow  rate.  The  relatively  small  number  of  data 
points  and  absence  of  reliable  information  near  the  wall,  contributed  to 
a  discrepancy  between  the  Re  calculated  from  the  velocity  profile  and 
measured  with  a  rotameter  of  about  This  estimate  was  obtained  using 
the  fact,  that  the  flow  rates  were  independent  of  tne  flow  regime,  and 
compairing  the  results  obtained  in  laminar  and  turbulent  flows  with  the 
constant  rotameter  reading.  The  accuracy  of  this  measurment  did  not  ap¬ 
preciably  changed  by  changing  the  flow  regime  in  the  pipe  from  laminar 


to  turbulent. 
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CHAPTER  3 


EXPERIMENTAL  RESULTS 


3.1  Mean  Flow:  Steady  vs.  Pulsating  Velocities  and  Pressure 


£ 

Fig. 3.1  represents  tine  mean  velocity  profile  measured  by  an 
x-wire  probe  in  laminar  pulsating  flow  and  normalized  by  the  center-line 
velocity.  The  velocity  profile  presented  was  taken  at  Re=4000,  period 
T=  1 . 34  sec,  relative  amplitude  of  bulk  velocity  U.j/U^O*.  In  Fig.3.1b 
three  measured  turbulent  velooity  profiles  are  shown,  one  of  the  pro¬ 
files  was  taken  in  steady  flow  at  Re*4000,  while  for  the  other  two  the 
amplitude  of  pulsations  was  varied.  No  significant  difference  can  be 
observed  between  steady  and  pulsating  time  mean  velocity  profiles  wheth¬ 
er  in  laminar  or  turbulent  flow.  This  result  could  be  anticipated  in 
fully  developed  laminar  flow,  beoause  of  the  resulting  linearity  of  the 
Navier-Stokes  equations,  but  in  turbulent  pulsating  flow  it  indicates, 
that  the  time  mean  Reynolds  stresses  are  not  affected  by  the  oscilla¬ 
tions  . 


The  frlotion  ooeffiolent  A  ,  calculated  from  Daroy's  formula 
P«  Al/D  U*/2 
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did  not  show  any  difference  between  steady  and  pulsating  flows  for  ootn 
laminar  and  turbulent  regimes.  The  measured  friction  coefficient  in 
laminar  flow  was  higher  by  aoout  1Q»,  than  the  tneoretically  predicted 
value  of  64/he.  In  turbulent  flow  the  value  of  \  obtained  is  in  fair 
agreement  with  the  values  quoted  in  the  literature  for  smooth  pipes  (see 
n.ochlichting  (1975)). 

A  number  of  factors  might  have  influenced  the  accuracy  at  wnien 
yjwas  defined.  1)  The  accuracy  of  the  measuraent  of  the  mean  bulk  velo¬ 
city,  as  discussed  above.  2)Pressure  was  measured  in  the  entrance  noz¬ 
zle  rather  than  in  the  developed  region  of  the  pipe  in  order  to  increase 
the  resolution  of  the  measuraent.  Pressure  differences  at  He=4000  along 
tne  entire  pipe  were  approximately  2  mm  of  water  in  laminar  flow  *.na  k 
mm  of  water  in  turbulent  flow.  The  measured  pressure  therefore  included 

the  influence  of  the  developing  flow  in  the  entrance  region  and  tne  dy- 

—2 

namic  head  component  pU  /2.  5)  It  is  difficult  to  null  the  output  of 
the  pressure  transducer  because  of  relatively  low  values  of  mean  pres¬ 
sure,  as  compared  to  the  pick  pulsation  values. 


3.2  Phase  Mean  Values:  Laminar  vs.  Turbulent  Plow 


Phase-averaged  data  is  obtained  from  the  measuraents  discussed  for 
ootn  laminar  and  turbulent  flow  regimes  thus  providing  the  first  two 
terms  in  the  decomposition  (1.1).  flg.3>2  shows  a  typical  dependence  of 
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phase  averaged  velocities  in  a  tuoulent  pulsating  flow  during  one  peri¬ 
od.  Tne  concomitant  reference  pressure  oscillations  are  snown  at  tne 
top  of  tne  figure,  tacn  velocity  trace  is  normalized  on  tne  time  .dear, 
velocity  on  tne  center  line  of  the  pipe,  ine  uppermost  velocity  trace 
corresponds  to  r/R=0  (i.e.  tne  velocity  was  measured  on  tne  center 
line)  and  the  oottom  trace  represents  tne  velocity  at  r/n=0.?7.  during 
a  fraction  of  tne  period  the  pressure  at  tne  inlet  of  tne  pipe  is  lower 
than  at  the  exit,  pointing  to  tne  existence  of  an  adverse  pressure  gra¬ 
dient.  The  velocity  however  does  not  reverse  itself  at  all  radial  posi¬ 
tions  and  at  all  pna3e  angles.  The  relations  oetween  the  amplitudes  an-, 
the  phase  angles  of  pressure  and  velocity  oscillations  nay  thus  lead  to 
i  situation  wnen  at  a  portion  of  tne  period  tne  direction  of  the  flow  is 
opposite  to  the  Instantaneous  direction  of  tne  pressure  gradient. 


The  validity  of  the  assumption  tnat  pressure  and  vel'bcity  oscilla¬ 
tions  are  harmonic  may  he  checked  6y  representing  tne  ensemore,  averaged 

\ 

signals  in  Fourier  series.  Tne  "power"  spectra  were  calculated,  ahp  tne 

\ 

ratio  of  the  two  first  coefficients  c(2  u>  )/c(  u> ) ,  i.e.  for  the  coeffi¬ 
cient  of  the  fundamental  frequency  and  its  first  harmonic,  was  deter¬ 
mined.  for  moderate  amplitudes  of  velocity  pulsations,  this  ratio  was 
less  than  3fr  thus  providing  the  justification  for  tne  harmonic  assump¬ 
tion.  At  high  amplitudes,  In  turbulent  pulsating  flow,  relaalnarlzation 
occurs  making  the  contribution  of  the  second  harmonic  more  significant. 

3.2.1  Oscillating  Part  of  the  Velocity  Profile 
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Phase  locked  velocity  profiles  can  Oe  obtained  from  the  phase  aver¬ 
aged  velocity  records.  Since  the  time  mean  velocity  profiles  are  inde¬ 
pendent  of  the  pulsations,  it  is  more  instructive  to  plot  the  oscillat¬ 
ing  part  of  the  velocity  profile  only  (i.e.  the  part  represented  by  the 
second  term  in  equation  (1.1)).  In  the  laminar  case  a  comparison  with 
tne  theoretical  predictions  of  Uchida  is  performed.  Although  at  least 
1024  points,  sampled  in  each  period,  made  it  possible  to  plot  phase 
locked  velocity  profiles  at  0.3°  intervals,  it  suffices  to  plot  eight 
profiles,  with  phase  difference  of  45°  for  the  purpose  of  the  following 
discussion. 

In  laminar  flow  the  shape  of  the  velocity  profiles  depends  only  on 
frequency  parameter  a  .  The  oscillating  part  of  the  laminar  velocity 
profiles  corrsponding  to  pulsation  periods  TsO.76  sec  (a  =11. d)  and 
T=2. 4  sec  (  a=6.7),  respectively,  are  shown  in  Fig.3.3a  and  Fig.3.3b. 
The  symbols  represent  the  measured  values  of  the  oscillating  velocity 
while  the  solid  lines  represent  the  theoretically  calculated  profiles, 
normalized  and  matched  to  the  measured  velocity  on  the  center  line.  The 
agreement  with  the  theory  is  good.  It  may  be  seen,  that  for  short  peri¬ 
ods  of  pulsation  the  flow  is  quite  uniform  in  the  oentral  core  of  the 
pipe  at  all  phases  of  the  cycle.  Sharp  velocity  gradients  occur  near 
the  wall.  The  extent  of  the  wall  region  is  proportional  to  />/<!).  This 
region  is  referred  to  in  the  literature  as  the  Stokes  layer  in  view  of 
the  analogy  with  the  oscillating  plane  boundary  layer  whioh  was  analysed 
by  Stokes.  For  large  a  the  Stokes  layer  beocmes  narrow  causing  the  vel¬ 
ocity  gradients  to  lnorease. 
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In  the  turbulent  flow  important  parameter  is  the  ratio  between  the 
thickness  <5  of  the  Stokes  layer  and  viscous  suolayer  6  .  If  one  as- 

sumes  that  the  effects  of  pulsations  are  limited  to  the  Stokes  layer 
while  the  effects  of  turbulence  are  excluded  from  the  viscous  sublayer 
(i.e.  from  0<yu./v  <5)  then  for  £  < 1  tile  pulsating  part  of  the 

O  V 

velocity  profile  in  laminar  and  turbulent  flows  snoulo  De  identical. 
This  hypothesis  may  be  checked  by  assuming  that  the  thickness  of  the 
viscous  subiayer  is  equal  to  -5  v/u# ,  and  tne  friction  velocity  u# 
may  for  the  sake  of  convenience  be  related  to  the  mean  velocity  U  by  tne 
power  law  (see  Schlichting  p.50tt) 

u,2=0.0225U7/4(v  /R)1/4 

one  obtains  the  following  expression  for  ''v  : 

6V  =5  v  /u#=60R/Re7/y 

which  renders  the  desired  thickness-ratio  between  of  tne  two  layers 

6  /  <5  =He7/y/60  a 

O  -  V 


At  Re=4000  6^/R=0.04  thus  for  T=0.7&  sec,  6  /  =2.0,  and  the 
effects  of  oscillations  on  the  radial  velocity  distribution  are  visible 
(Pig. 3. 4  ).  At  lower  frequencies  turbulence  penetrates  into  the  Stokes 
layer  and  destroyes  it,  and  the  frequency  parameter  a  is  no  longer  im- 
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portant  resulting  in  a  much  more  uniform  velocity  profile  (Pig. 3.4  re¬ 
corded  at  '1=2.4  sec,  giving  6 ^ /6y  =3.5.  In  order  to  reduce  furtner 
the  thickness  ratio  6  /6^  tne  ne  was  reduced  to  2900  for  T=0.55  sec 
resulting  in  6^/6^  =1.2.  Tne  measured  oscillatory  velocity  profiles 
arc  plotted  in  Pig. 3. 4°.  There  is  a  reasonaDie  similarity  between  the 
measured  turbulent  (symbols)  and  calculated  theoretically  laminar  (solid 
line)  oscillating  velocity  profiles  for  this  case. 

3 .2.2  fladial  Distribution  of  Velocity  Phase  Angle  and  Amplitude 


The  harmonic  character  of  pulsations  enables  one  to  represent  the 
instanteneous  velocity  by  an  exponential  form  (tiq.( 1 .3)) ,  and  aleviates 
the  necessity  of  describing  the  temporal  and  spatial  changes  in  velocity 
by  plotting  a  large  number  of  pnase  averaged  velocity  profiles.  Two 
functions  can  fully  describe  the  oscillating  component  of  velocity  at 
tne  imposed  frequency:  (i)  the  amplitude  distribution  u^r);  (li)  the 
phase  angle  $u(r)  relative  to  the  phase  of  the  pressure  oscillations. 
In  order  to  obtain  these  functions,  tne  phase  averaged  velocity  was 
Fourier  transformed 


* 

u(r,t)»a  (r)+yU  (r)cos(nt/T)+b(r)sin(nt/T))  3-1 

o  *•.  n  n 

n»l 


where  N  is  half  of  the  number  of  points  sampled  at  each  period,  and  T  is 
period  of  pulsations.  The  coefficients,  that  describe  the  oscillations 
at  the  forcing  frequency  1/T,  were  then  used  to  find  the  amplitude  fuc- 
tion  u^rJs/a^rJ+b^Cr)  and 


the 


initial 


phase 


angle 


and 
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♦  (r)=arctan((b,tr)/a,(r)) .  from  tne  initial  pnase  of  the  velocity  os- 
o  II 

dilations,  tne  initial  pnase  of  tne  pressure  oscillations  was  suo- 
stracted,  giving  phase  angles  $u(r)  relative  to  tne  pressure. 

in  tne  aoser.ce  of  viscosity,  tne  pressure,  oeing  tne  only  driving 
force,  is  in  phase  with  the  acceleration  of  the  fluid.  Ihe  velocity 
lags  oy  90°  behind  the  acceleration,  and  thus  also  lags  oenind  tne  pres¬ 
sure.  The  radial  distribution  of  the  phase  angle  is  shown  in  Fig. 3.5 
for  fie»4000  and  various  periods  of  oscillations.  The  solid  lines  show 
the  theoretical  prediction  of  Uchida,  while  the  crosses  and  the  trian¬ 
gles  give  the  measured  pnase  angles  In  fully  developed  laminar  and  tur¬ 
bulent  flows,  respectively,  k  good  agreement  with  tne  theory  was  obta¬ 
ined  in  the  laminar  case,  with  tne  exception  of  very  low  frequencies  for 
which  the  influence  of  the  entrance  region  becomes  more  pronounced,  as 
it  will  be  described  later.  The  phase  lag  on  the  center  line  in  laminar 
flow  is  usually  90°,  and  it  decreases  to  approximately  45°  near  the 
wall . 


In  turbulent  flow  the  results  appear  to  be  very  different.  The 
frequency  parameter  a  *R/ui/v  no  longer  controls  the  flow  because  the  re¬ 
levant  viscosity  becomes  some  Kind  of  turbulent  exchange  coefficient  e  , 
which  is  orders  of  magnitude  larger  than  v  .  The  effective  frequency 
parameter  a  is  thus  much  lower,  and  the  phase  lag  of  the  velocity  in  the 
central  region  of  the  pipe  decreases  more  quickly  with  Increasing  the 
period  than  in  the  corresponding  laminar  flow.  In  contrast  to  the  lami¬ 
nar  flow,  the  phase  lag  increases  towards  the  wall.  The  qualitative  na- 
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ture  of  the  result  .-as  noticed  oy  namaprian  and  tu  i  I .  «  more  deta¬ 
iled  explanation  oi  tnis  pneno.oer.or*  wnicn  is  o\suj  on  a  simple  turoulent 
model  for  the  pulsating  pipe  flow  will  oe  discusses.  later. 

The  radial  distribution  of  the  phase  lag  in  velocity  in  the  turbu¬ 
lent  flow  is  both  amplitude  and  mean  He  dependent.  The  dependence  on 
the  amplitude  of  pulsations  is  rather  weax.  The  radial  distribution  of 
the  phase  lag  for  two  amplitudes  of  pulsations  at  two  frequencies  is 

ci  b 

presented  in  figures  3.6  and  3.0  .  For  longer  periods  of  pulsations 
(T=2.4  sec)  changes  in  amplitude  cause  no  concomitant  change  in  the 
central  region  of  the  pipe.  At  higher  frequencies  (T=1.25  sec)  the  in¬ 
fluence  of  the  amplitude  on  the  phase  angle  is  felt  accross  the  entire 

pipe,  hevertneless ,  a  large  increase  in  the  amplitude  causes  5  rela- 

o 

tively  small  change  in  the  phase  angle  (about  .1  ). 

Tne  dependence  of  the  phase  angle  on  the  mean  heynolds  number  is 
more  pronounced.  A  representative  viscosity  increases  with  increasing 
Re,  reducing  effective  a,  Tne  velocity  phase-lag  therefore  also  decre¬ 
ases  with  increasing  Re  (fig. 3.7).  On  the  other  hand,  the  viscous  sub¬ 
layer  becomes  thinner  with  increasing  Re,  and  the  radial  distribution  of 
the  phase-lag  at  Re=7500  is  practically  constant. 

A  qualitative  difference  between  laminar  and  turbulent  pulsating 
flows  also  exists  In  the  radial  distribution  of  the  amplitudes  of  the 
velocity  oscillations.  In  laminar  flow,  the  maximum  amplitude  of  the 
velocity  oscillations  occurs  In  the  Stokes  layer  near  the  wall,  as  no- 


t 
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ticed  by  Ricnardson  (1cj2f/2o),  v.hile  in  turbulent  flow  tne  maximum  am¬ 
plitude  occurs  in  the  center  of  tne  pipe,  lne  radial  distribution  of 
the  measured  amplitudes  of  velocity  oscillations  in  laminar  and  turbu¬ 
lent  flows  resulting  from  identical  forcing  is  snown  in  Fig. 3 .d.  The 
theoretically  calculated  amplitudes  in  laminar  flow  which  were  matched 
on  the  center  line  (Uchida  (1^56))  are  shown  also  in  Fig. 3. 6  for  compar¬ 
ison. 


The  amplitude  distribution  of  the  axial  component  of  velocity  in 
the  turbulent  flow  is  nearly  uniform  in  the  central  region  of  the  pipe, 
out  decreases  rapidly  near  the  wall.  The  amplitude  of  pulsations  in 
laminar  flow  also  decreases  near  the  wall,  increasing  the  mean  He  or 
decreasing  the  frequency  of  oscillations  lead  in  turbulent  flow  to  more 
uniform  distribution  of  the  velocity  amplitudes  for  the  reasons  already 
discussed  in  conjunction  with  the  radial  distribution  of  the  phase  lag. 

3.2.3  The  Influence  of  the  Entrance  Region 

The  disagreement  between  the  calculated  and  measured  pnase  angles 
of  velocity  at  low  frequencies  may  be  attributed  to  the  influence  of  the 
entrance  region  on  the  flow.  To  check  this  further,  a  second  pressure 
transducer  was  placed  100  diameters  upstream  the  exit  of  the  pipe,  (i.e. 
400  diameters  from  the  first  one,  which  was  located  at  the  entrance  of 
the  pipe).  The  flow  at  x/Ds400  is  considered  to  be  fully  developed, 
thus  tne  comparison  between  the  two  measured  pressures  could  provide  in¬ 
formation  about  tne  development  of  pressure  gradient  along  the  pipe. 
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phase  shift  in  prESSuro  pulsations  was  found  to  exist  between 
tneu-  two  positions  and  is  shown  in  fig.o.S11  as  a  function  of  the  pulsa- 
tion  period  in  botn  laminar  and  turoulent  flows,  ho  notable  amplitude 
dependence  was  found  in  the  aeasurments  whian  were  made  at  mean  Re=4000. 
«.t  long  periods  of  forcing,  cressure  at  the  exit  leads  tne  pressure  pul¬ 
sations  at  tne  inlet  of  the  pipe.  The  phase  difference  is  much  nigher 
in  iaminar,  than  in  turbulent  flow,  and  decreases  in  both  cases  with  de¬ 
creasing  period  of  pulsations. 

The  calculated  values  of  tne  same  phase  shift  are  3hown  also  in 
t  lgo.9a.  The  calculations  were  made  for  the  appropriate  values  of  Re 
and  T  according  to  the  theoretical  calculations  of  Atabek  and  Cnang 
(1*61)  for  the  laminar  pulsating  flow  in  the  entrance  region  of  a  pipe, 
it  is  clearly  seen  that  the  phase  shift  angles  are  underestimated  by  tne 
theory.  These  results  support  the  conclusion  of  Denison  (1970),  wno 
found  that  the  theory  of  Atabek  and  Chang  predicts  a  faster  evolution  of 
the  pulsating  flow  to  the  asymptotic,  fully  developed  form,  than  ob¬ 
served  experimentally.  The  conclusion  of  Denison  is  based  on  measur- 
ments  of  pulsating  velocity  profiles  in  the  iaminar  entrance  region  of 
tne  pipe.  It  should  be  noted,  however,  that  qualitatively  the  predic¬ 
tion  is  correct:  the  pressure  at  inlet  lags  behind  the  pressure  in  the 
fully  developed  region,  and  the  pnase  difference  increases  with  increas¬ 
ing  the  period  of  pulsations.  The  influence  of  the  entrance  region  in 
tne  turbulent  pulsating  flow  is  much  weaker,  than  in  laminar  flow,  and 
phase  differences  are  therefore  less  pronounoed. 
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At  higher  frequencies  (T<0.7  sec),  the  sign  of  the  pnase  difference 
cnanges.  Tne  finite  value  of  the  sound  velocity  becomes  important  at 
these  frequencies,  and  the  pressure  downstream  lags  behind  the  cnanges 
in  upstream  pressure  as  predicted  in  Section  2.5.  The  measured  phase 
differencies  are  slightly  higher  than  the  theoretically  calculated.  The 
discrepancies  may  stem  from  the  finite  angular  resolution  in  the  experi¬ 
mental  setup  (about  0.3°  for  each  transducer)  and  the  hydraulic  approxi¬ 
mation,  used  in  the  theory.  It  is  interesting  to  notice,  that  there  is 
practically  no  difference  in  the  measured  pressure  phase  lags  between 
laminar  and  turbulent  pulsating  flows  at  tnis  range  of  frequencies  (T<1 
sec) . 


An  additional  way  to  check  the  influence  of  the  entrance  region  is 
to  compare  the  amplitudes  of  the  pressure  pulsations  at  different  stre- 
amwise  locations.  As  it  was  mentioned  earlier,  one  transducer  was 
placed  at  inlet  of  the  pipe  which  is  500  D  long,  while  the  second  trans¬ 
ducer  was  placed  100  D  from  the  exit.  Thus  for  a  linear  distribution  of 
pressure  along  the  pipe  the  ratio  of  the  pressure  amplitudes  should 
therefore  be  equal  to  0.2.  This  ratio  would  be  reduced  as  a  result  of 
the  influence  of  the  entrance  region.  The  accuracy  of  measuring  the  am¬ 
plitude  is  much  better  than  the  accuracy  in  measuring  the  absolute  stea¬ 
dy  pressure  because  the  former  is  independent  of  the  error  in  the  zero 
setting. 

Flg.3.9b  showes  the  measured  ratio  of  the  pressure  amplitudes  in 
both  laminar  and  turbulent  flows.  There  is  no  significant  difference 
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oetween  these  flow  regimes;  the  results  suggest,  however,  that  the  in¬ 
fluence  of  tno  entrance  region  diminishes  with  increasing  frequency.  It 
is  worth  noting  that  the  effects  of  the  entrance  region  are  much 
stronger  on  the  steady,  than  on  the  oscillating  parts  of  the  flow.  This 
stems  from  the  fact  that  the  oscillating  velocity  profile  in  fully  de¬ 
veloped  laminar,  or  turbulent  flows,  does  not  differ  significantly  from 
the  slug-type  velocity  profile,  which  enters  the  pipe  fran  the  settling 
chamber.  The  theory  of  Atabek  and  Chang  does  not  predict  notable  nonli¬ 
nearity  in  the  distribution  of  the  pressure  amplitude  along  the  pipe. 

3.2.4  The  Relation  between  the  Pulsations  of  Pressure  and  Flow  Rate 

The  mechanical  power  U,  neccessary  to  push  the  flow  at  a  rate  Q 
through  a  pipe  in  which  the  pressure  drops  by  Ap  is  proportional  to  the 
product  Q  Ap.  In  pulsating  flow  the  power  is  time  dependent,  and  the 
amplitude  of  power  oscillations  is  determined  by  the  amplitudes  and 
phase  relations  between  the  pulsations  of  the  pressure  and  flow  rate: 
WspjQ^cos  ^ .  In  order  to  compare  the  amplitudes  of  power  pulsations  in 
laminar  and  turbulent  flows  the  relation  between  p 1  and  Q1  has  to  be 
known. 

Fig. 3. 10  shows  the  dependence  of  the  measured  amplitude  of  pulsa¬ 
tions  of  bulk  velocity  on  the  amplitude  of  the  imposed  pressure  oscilla¬ 
tions  for  two  frequencies.  Only  small  differences  can  be  noticed  in  the 
responce  of  the  velocity  amplitude  to  the  imposed  forcing  in  laminar  and 
turbulent  flows. The  bulk  velocity  amplitude  is  proportional  to  that  of 
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pressure,  and  tne  difference  Between  the  slope  of  the  laminar  and  the 
turbulent  dependence  is  within  the  experimental  error.  For  a  given  am¬ 
plitude  of  pressure  pulsations  tne  amplitude  of  the  bulk  velocity  is 
nigner  at  lower  frequency  as  it  follows  from  the  expected  dependence  of 
<r  on  the  frequency  (see  Fig. 2. 2). 

q 

At  a  given  frequency  of  pulsations  only  the  phase  difference  4> 
can  contribute  to  the  difference  in  the  amplitude  of  power  pulsations 
between  the  laminar  and  the  turbulent  flows .  In  laminar  flow  tne  phase 
angle  does  not  differ  appreciably  from  90°  (it  changed  from  about  b5°  at 
T=0.56  sec  to  about  70°  at  T=2.4  sec).  In  the  turbulent  flow,  changes 
in  frequency  cause  a  more  significant  decrease  in  the  phase  lag  of  the 
flow  rate  behind  tne  pressure,  as  it  qualitatively  can  be  concluded  from 
Fig. 3.5.  At  high  frequencies,  however,  there  is  only  a  small  difference 
in  ^  between  laminar  and  turbulent  flows;  at  1=0.56  sec  it  is  about 
05°  in  both  cases.  In  turbulent  flow  decreases  fast  with  decreasing 
frequency  (at  T=2.4  sec  <^=50°),  and  cos  therefore  becomes  larger, 
than  in  the  corresponding  laminar  flow. 

There  is  thus  no  significant  difference  in  the  amplitude  of  the 
power  pulsations  between  laminar  and  turbulent  flows  at  high  frequen¬ 
cies,  where  the  oscillating  velocity  prifiles  are  very  much  alike  (see 
Section  3.2.1).  Decreasing  the  frequency  increases  the  amplitude  of  the 
flow  rate  for  given  pulsations  of  pressure  and  the  amplitude  of  the 
power  pulsations  grows  in  both  laminar  and  turbulent  flows;  the  growth 
in  the  turbulent  case  is  faster  because  of  stronger  dependence  of  on 
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frequency. 


3-3  Turoulent  Characteristics 


There  is  no  doubt,  that  a  lot  of  information  about  the  structure  of 
the  turbulent  flow  is  lost  by  the  conventional  and  phase  locked  averag¬ 
ing  procedure.  In  Fig. 3. 11  the  instantaneous  axial  u  and  radial  v  velo¬ 
cities,  measured  during  a  single  period  are  shown  and  compared  with  the 
phase  averaged  values;  the  dependence  of  the  pressure  on  time  is  shown 
in  the  upper  curve.  Several  conclusions  could  be  drawn  from  this 
Figure:  1)  the  phase  averaged  value  of  the  radial  velocity  vanishes,  so 
there  is  no  Reynolds  stress  component,  connected  with  the  orderly  pul¬ 
sating  part  of  flow  as  expected  in  the  fully  developed  pipe  flow;  2) 
the  intensity  of  radial  velocity  fluctuations  is  of  the  same  order  of 
magnitude,  as  the  intensity  of  axial  velocity  fluctuations;  3)  there  is 
some  evident  dependence  of  the  turbulent  activity  on  the  phase  of  im¬ 
posed  oscillations  which  is  expressed  in  the  changes  in  both  amplitude 
and  frequency  of  the  turbulent  velocity  fluctuations  (see  also 
Fig. 3. 12). 

The  turbulent  activity  becomes  increasingly  dependent  on  the  phase 
of  forcing  at  higher  amplitudes  of  the  pulsations.  The  oscillating  part 
of  velocity  during  a  single  period,  measured  by  the  hot-wire  rake  and 
compared  with  the  signature  of  the  pressure  oso illations,  is  shown  in 
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fig. j. 12.  A  partial  iaminarization  of  the  flow  is  ooserveci  during  tne 
tine  corresponding  to  tne  minimum  velocities  produced  by  tne  forcing,  k 
" oreaKcio’/.n"  to  occur  at  a  tine  corresponding  to  tne  maximum  in- 
stanteneous  velocity,  wnen  tne  aaplitude  of  turoulcnt  fluctuations  in¬ 
creases  quite  sudcenly.  The  decrease  in  tne  amplitude  of  tne  fluctua¬ 
tions  on  the  decelerating  portion  of  the  cycle  is  relatively  gradual. 

At  even  higher  amplitudes  of  pulsations  the  flow  becomes  completely 
laminar.  At  these  amplitudes  the  instantaneous  Re  falls  below  its  crit¬ 
ical  value  during  a  fraction  of  the  period.  Vrnen  the  duration  of  tne 
flow  at  tne  suocrltical  Re  is  sufficiently  long,  a  complete  relaminari- 
zatlon  occurs.  The  mean  velocity  of  the  flow  at  Re=4000  is  about  l.o 
m/sec,  so  that  even  at  lowest  frequency  of  pulsations  (T=4.p  sec)  each 
fluid  particle  remained  in  tne  pipe  for  several  periods.  The  flow  is 
naturally  laminar  at  ne=4000,  and  tne  transition  to  turbulence  was  trig¬ 
gered  artificially  oy  a  perturbance,  which  was  purposely  placed  at  the 
entrance  to  the  pipe.  Thus  if  the  flow  relaminarizes  further  downstre¬ 
am,  it  will  remain  laminar  throughout  the  pipe.  This  effect  was  noticed 
by  a  number  of  investigators,  cited  earlier  (e.g.  Sarpxaya  (1966)  and 
hamaprian  and  Tu  (I90O)). 

3.3*1  The  Intensity  of  Velocity  fluctuations 

The  distribution  of  the  time  averaged,  root-mean  squared,  turbulent 
fluctuations  in  the  axial  and  radial  directions  for  steady  and  pulsating 
flows  is  compared  in  fig. 3. 13.  The  velocity  fluctuations  are  normalized 
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oy  the  friction  velocity  uu  wnicn  was  calculated  from  tne  time  mean 
value  of  the  pressure  drop  (u^slp.l  cm/ sec  at  tte=4000).  There  is  no 


.2 


JP: 


significant  difference  oetween  time  mean  values  of/u*  and/v'  in  steady 

? 

and  pulsating  flows,  A  small  Out  conslstant  decrease  in  v'  can  be  no- 


o 

ticed  (fig. 3- 13  ),  with  increasing  amplitudes  of  pulsations.  The  dis- 


and 


trioution  of/u'*"  andPv,w  is  in  good  agreement  with  the  results,  obtained 
for  steady  turbulent  pipe  flow  by  Laufer  ( 195*1)  and  by  wygnanski  and 
Cnampagne  (1975). 


measurments  of  Laufer  were  made  at  AespOOGO,  which  is  higher  by 
order  of  magnitude,  than  the  prevailing  Re  used  in  the  present  investi¬ 
gation.  Nevertheless,  the  comparison  of  the  normalized  results  is  still 
possible,  because  Rho  values  of  the  turbulent  velocity  fluctuations  nor¬ 
malized  by  tne  friction  velocity  are  practically  independent  of  Re.  On 
the  other  hand,  uw  is  dependent  on  He,  and  it3  in3tanteneous  values 
change  in  the  pulsating  pipe  flow.  Thus,  u'  and  v*  are  pnase  dependent. 
This  conclusion  was  reached  by  examining  the  recorded  velocity  depen¬ 
dence  on  time,  as  shown  on  figures  3.11  and  3.12. 


The  phase  dependence  of  the  RMS  values  of  velocity  fluctuations  is 
plotted  in  Fig. 3. 14,  after  tne  calculated  data  was  smoothed  slightly  by 
a  running  average  procedure.  The  smoothing  was  necessary  because  of  re¬ 
latively  short  length  of  the  ensamble  used  in  the  averaging  process  (200 
measured  data  points  at  each  phase  angle).  A  number  of  conclusions  can 
be  drawn  from  the  examination  of  this  Figure:  1)  the  RMS  values  of 
axial  and  radial  velocity  fluctuations  are  pnase  dependent;  2)  the 
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phase  angle  of  the  Intensity  of  the  axial  and  radial  velocity  fluctua- 

2  2 

tions  is  not  identical;  3)  tne  pnase  dependence  of  /<ir  >  and  ^  <v%e”> 
is  obviously  non-harmonic ,  the  growth  is  notably  faster  than  the  decay, 
as  it  was  noted  previously  wnile  examining  the  instanteneous  velocity 
records  over  a  single  period. 

— T 

In  order  to  compare  the  radial  dependence  of  the  phase  of  /<u’p 
and  \j  <v’^>  with  tne  phase  of  the  mean  velocity,  it  was  decided  to  limit 
the  analysis  to  the  fundamental  frequency  only,  inspite  of  tne 
non-harmonic  behaviour  of  this  ensemble  averaged  data. 

2~ 

Figure  3.15  shows  the  phase  lag  of  /  <u'  >  behind  the  pressure  and 
compares  it  to  the  phase  of  the  mean  velocity,  for  three  different  fre¬ 
quencies  of  pulsations  Dut  approximately  for  identical  amplitude  of  the 

bulk  velocity.  The  higher  the  frequency,  the  larger  is  the  phase  lag  of 
- T 

/  <u'  >  behind  tne  fundamental  oscillation  of  the  flow  velocity.  It  is 

2~ 

interesting  to  note,  that  the  phase  angle  of  /<u’  >  changes  signifi¬ 
cantly  with  the  change  in  the  radial  position,  being  maximum  at  the 
center  of  the  pipe,  and  attaining  a  minimum  in  the  region  corresponding 
to  the  maximum  turbulent  production.  At  lower  frequencies  the  pulsa¬ 
tions  (T>2  sec)  /  <u'2>  at  r/fl=0.7  may  even  lead  the  phase  averaged  vel¬ 
ocity.  In  Fig. 3.1^  the  phase  distributions  of  <u>  and  /<u'2>  are  com¬ 
pared  for  a  constant  frequency  of  forcing  (1=2.4  sec)  and  the  Reynolds 
numoers  ranging  from  3300  to  7500.  The  pnase  difference  oetween  <u>  and 

- 7 

/ <u '  >  decreases  with  increasing  he;  the  location  of  minimum  in  the 
phase  difference  approaches  the  wall  at  high  values  of  the  Reynolds 
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number . 

The  phase  angle  of  /<v,<L>  is  practically  independent  of  tne  radial 
position  (Fig.i.lo).  The  phase-lag  of  /<v'^>  behind  pressure  is  very 
close  to  that  of  /<u’ 2>  in  tne  center  of  tne  pipe ,  Dut  it  does  not  de¬ 
crease  at  larger  r/h.  The  magnitude  of  tne  phase-lag  in  /<u,2>  and 
— 2" 

Av'  >  relative  to  the  pulsating  velocity  depends  on  the  amplitude  of 
the  velocity  pulsations,  but  the  radial  distribution  of  the  phase  is 
unaffected . 

— 2“ 

Tne  qualitative  difference  in  the  phase  benaviour  of  /<u'  >  and 

2~ 

Av'  >  can  be  explained  by  considering  the  oudget  of  the  turbulent  en¬ 
ergy  (see  hinze,  p.325).  The  rate  of  change  in  the  intensity  of  the 

2 

longitudinal  fluctuations  (i.e.  3<u'  >/  3 1 )  is  governed  by  the  produc¬ 
tion  term  -<u'v’>  3<u>/  3r,  so  that  simultaneous  existence  of  the  oscil¬ 
lating  velocity  gradient  and  Reynolds  stress  is  necessary  for  the  pro¬ 
duction,  as  it  was  pointed  out  by  Rotta  (1962).  The  temporal  change  in 
the  phase-mean  velocity  gradient  produces  a  most  pronounced  change  in 
✓<u*  >  at  the  location  where  the  product  of  the  oscillating  Reynolds 
stress  with  3<u>/  3r  attains  maximum.  The  phase  difference  between  <u> 

"  p-" 

and  /<u'  >  is  thus  minimal  in  the  region  of  maximum  production  and  oc¬ 
curs  approximately  at  r/Rs0.7  at  time-mean  Res4000;  at  higher  values  of 
Re  the  location  of  the  minimum  in  the  phase  difference  approaches  the 
wall.  Once  generated,  the  longitudinal  velocity  fluctuations  are  trans¬ 
ferred  acoross  the  pipe  by  energy  redistribution  mechanisms,  which  are 
connected  to  pressure  fluctuations.  The  phase  lag  between  <u>  and 
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ii:cr-;  5e„  therefor,  wit..  t..a  increasi..;,  ic.tar.cc  from  tna  region 
iia.:  proa  motion.  ;.i.. :.c  ^ licet  pr  >«.  -atior  j;  tne  turouler.t 
...'j;';;  ir,  tr.c.  central  re; ion  oi  tr.c  pi,-.-  resultin'  in  i  mcjcimum  pause 
i  .  taere . 

fhe  radial  velocity  fluctuations  do  not  extract  energy  directly 

2 

from  the  mean  flow  but  rather  from  <u*  >  through  pressure  redistribution 
terms.  Thus  the  phase-lag  in  /  <v,2>  is  similar  to  the  phase-lag  of 
/<u,2>  in  the  center  of  tne  pipe,  out  the  pnase-lag  of  /  <v'^>  is  inde¬ 
pendent  of  radial  location. 

in  addition  to  the  pnase  oehaviour  of  tne  «M£>  of  tne  turbulent  vel¬ 
ocity  fluctuations,  the  racial  distrioution  of  tneir  amplitudes  was 
onoc.ccd  as  well,  ihe  dependencies  of  tne  dimensionless  amplitudes  of 

p 

Au'  >,  normalized  by  the  amplitudes  of  tne  pulsating  ouIk  velocity  for 
tnree  periods  of  forcing,  are  plotted  in  fig. 5. 17.  increasing  the  fre¬ 
quency  of  the  pulsations  results  in  higher  dimensionless  amplitudes  of 
/<u'S.  A  large  change  in  the  amplitude  of  tne  bulx  velocity  produces 
a  small  growth  in  the  amplitude  of  /<u'  >.  The  data  presented  in 
fig. j. 1b  reflects  a  change  of  approximately  30Q>  in  tne  amplitude  of  tne 
ouIk  velocity  pulsations  which  lead  to  a  cnange  of  20}  in  tne  dimension¬ 
less  amplitude  of  /<u,2>. 

in  order  to  check  the  non-harmonic  benaviour  of  /<u'^>  and  /<v,2> 
tne  ratio  of  tne  first  two  power  spectral  coefficients  c(2  o>)/c(  u>)  was 
calculated  for  tne  pnase  averaged  mm*  values.  The  first  coefficient 


Page  59 


c( ui )  corresponds  to  the  fundamental  frequency  of  pulsations,  while 
c(2<d)  to  its  first  harmonic.  No  significant  frequency  or  radial  depen¬ 
dence  of  c(2  on  )/c(  a) )  was  observed.  This  ratio  was  about  18>  for 
Uj/0=15J,  and  increased  to  c(2  to  )/c(  to  )=40>  at  U 1  /U=35% .  The 
non-harmonic  behaviour  of  /<v,Z>  was  similar.  The  phase  dependence  of 
v<u,Z>  is  becoming  more  harmonic  when  Re  is  increased.  At  Re=5800  and 
Re=7500  c(2 to  )/c( to )  is  less  than  5%  at  all  radial  locations.  The  con¬ 
clusion  thus  can  be  made  that  the  non-harmonic  behaviour  of  the  turbu¬ 
lent  characteristics  at  Re<H500  is  caused  by  the  onset  of  the  relaminar- 
ization  process  at  low  instanteneous  values  of  the  bulk  velocity. 


3.3.2  Reynolds  Stresses 


Time  averaged  Reynolds  stresses  -u'v'  which  were  measured  by  an 
x-wire,  are  essentially  independent  of  the  flow  pulsations  (Fig. 3. 19), 
as  is* the  case  of  the  RMS  values  of  turbulent  velocity  ^fluctuations. 
This  conclusion  was  inferred  earlier  from  the  similarity  of  time  mean 
velocity  profiles  (Section  3-1).  The  fact  that  the  time  mean  pressure 
drop  along  the  pipe  is  independent  of  forcing  also  indicates,  that  the 
pulsations  have  nearly  no  effect  on  the  mean  Reynolds  stresses. 


The  radial  distribution  of  the  oscillating  part  of  the  Reynolds 
stresses  for  two  amplitudes  of  pulsations  is  shown  in  Fig. 3. 20.  The 
measured  values  of  <u'v'>  were  rendered  dimensionless  when  divided  by 
the  amplitude  of  the  pulsations  of  uc  which  in  turn  were  deduced  from 
the  oscillatory  component  of  the  pressure  drop.  Normalized  in  that  way, 
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the  amplitude.*?  of  the  Reynolds  stress  appear  to  be  independent  of  the 
amplitude  of  the.  pulsating  bulx  velocity.  (qualitatively,  the  radial 
distribution  of  tne  amplitude  of  <u'v*>  resembles  the  time  mean  u’v'  de¬ 
pendence,  both  srow  linearly  with  tne  radius  in  the  central  region  of 
the  pipe,  but  tne  pulsating  <u*v'>  attains  its  maximum  value  closer  to 
the  wall,  than  the  time  mean  Reynolds  stresses. 

The  phase  lag  of  <u'v'>  behind  pressure  (Fig. 3.21)  is  similar  to 
the  phase-lag  of  Au,2>  at  comparable  amplitudes  and  radial  positions 
(compare  with  Fig. 3. 1b).  The  phase-lag  of  <u*v*>  is  minimal  not  far 
from  the  wall,  at  0.7<r/R<0.6;  with  the  location  of  minimum  depending 
on  the  amplitude  of  the  forcing.  The  reseoblence  between  the  radial 
distribution  of  the  phase-lag  of  /<u*  >  and  <u*v'>  stems  from  the  fact 
that  u'v'  fluctuations  can  extract  energy  directly  from  the  mean  flow 
through  the  production  term  -<v*  >  9<u>/  9r  (Hinze  p.325).  The  minimum 
in  the  radial  distribution  of  the  phase  angle  of  <u'v'>  is  much  less 
pronounced  than  for  the  longitudinal  velocity  fluctuations,  because  of 
the  uniform  distribution  of  both  amplitude  and  phase  angle  of^<v'  >. 

The  Reynolds  stresses,  as  the  RMS  values  of  velocity  fluctuations, 
are  not  changing  harmonically  in  time.  The  "power  speotrum"  of  the 
phase  averaged  <u*v'>  was  oaloulated,  and  the  coefficients,  correspond¬ 
ing  to  the  fundamental  frequenoy  of  pulsations  and  its  first  harmonic, 
were  compered.  Their  ratio  was  found  to  be  similar  to  that,  obtained 

— y  . 

for  ^u'  >,  l.e.  about  20%  when  the  amplitude  of  foroing  was  U^/UsZO} 
and  about  40»  for  Uj/Ua35l. 
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3-3-3  Spectral  Analysis 

One  of  tue  interesting  questions,  concerning  turoulent  pulsating 
pipe  flow,  is  the  interaction  of  turbulence  with  the  velocity  oscilla¬ 
tions  generated  by  forcing.  Cheng  (1971 )»  for  example,  states,  that 
there  is  a  strong  amplification  of  the  turbulent  fluctuations  at  the 
frequency  of  the  pulsations. 

In  order  to  check  the  spectral  distribution  of  the  energy  of  turbu¬ 
lent  fluctuations  at  frequencies  corresponding  to  the  frequency  of  the 
imposed  oscillations,  the  sampling  period  was  increased  to  eight  periods 
of  forcing.  The  data  was  acquired  in  the  following  way.  The  first 
trigger  pulse  recognized  by  the  computer  initiated  the  record.  From 
that  point  in  time  only  every  6-th  or  16-th  data  point  was  stored  after 
the  corresponding  phase-averaged  value  of  the  velocity  was  substracted 
from  it,  thus  leaving  only  the  fluctuating  turbulent  component  of  the 
instantaneous  velocity.  One  hundred  and  twenty  eight  (126)  data  points 
were  obtained  for  each  period,  giving  a  total  of  1024  points  per  record. 
The  power  spectra  were  calculated  for  each  record  and  averaged  over  the 
number  of  records  used. 

A  log-log  plot  of  the  speotra  at  9  radial  locations  is  shown  in 

Fig.  3  .22*  for  unforced  flow;  and  in  Fig. 3.22b  and  Fig. 3 .22°  for  foroing 

at  a  period  equal  to  T«1.34  seo  and  at  two  amplitudes  of  osoillations 

Uj/U*20i  and  1^/11*35 )>.  The  general  shape  of  the  speotra  of  the  steady 
and  pulsating  flows  is  very  similar.  Foroing  produces  a  strong  peak  in 
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tne  spectrum  at  the  frequency  of  tne  pulsations  in  tne  nonstationary 
flow,  in  tne  central  region  of  the  pipe  even  tne  second  and  tne  third 
narmonic  of  tne  pulsating  frequency  are  visible.  Tnis  confirms,  at  the 
first  signt,  tne  conclusion  of  Cheng  aoout  the  increase  of  tne  turbulent 
energy  at  the  frequency  of  the  pulsations,  To  check  the  matter  further 
the  same  type  of  measurment  was  repeated  in  laminar  pulsating  flow  at 
identical  flow  conditions;  tne  ensued  spectra  nad  identical  peaks  at 
the  frequency  of  the  pulsations  and  its  harmonics,  dinoe  the  intensity 
of  these  spectral  peaxs  relative  to  the  background  was  identical  to  the 
intensity  measured  In  turbulent  flow  and  yet  tne  flow  was  laminar,  error 
in  tne  processing  of  data  was  suspected. 

*s  explained  in  Section  2.5,  the  sampling  frequency  was  calculated 
oy  tne  computer  in  order,  that  1024  or  2046  data  points  would  be  sampled 
per  cacn  channel  per  period.  In  reality,  however,  the  desired  number  of 
sampled  data  points  per  period  was  not  precisely  achieved.  The  reason 
for  the  discrepancy  stems  from  the  fact  that  11  cnannels  were  recorded 
at  a  total  sampling  rate  of  approximately  15  Khz.  Inus,  about  70  usee 
elapsed  between  the  consecutive  data  points  wnile  the  resolution  of  tne 
cIock  being  only  1  usee,  resulted  In  a  maximum  possible  error  in  the 
sampling  frequency  of  approximately  1.5*.  The  sampled  period  might  have 
been  either  longer  or  shorter  then  the  real  period  of  pulsations,  re¬ 
sulting  in  a  disoontuinity  at  the  and  of  each  record.  The  absolute 
value  of  the  diaoontinulty  was  not  large,  and  it  oould  not  be  noticed 
from  observing  the  velocity  signature  on  a  "Tektronix"  screen,  but  its 
periodicity  caused  the  generation  of  the  spectral  peaks  discussed  in 
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conjunction  witn  fig. 3. 22  at  the  frequenoy  of  pulsations  and  its  harmon¬ 
ics. 


examination  of  the  low  frequency  spectra  excluding  the  above  menti¬ 
oned  3pikes  leads  to  the  conclusion,  that  there  is  no  significant  inter¬ 
action  between  the  imposed  oscillations  and  the  turbulence.  Thus  the 
conclusion  of  Cheng  might  be  in  error. 

In  order  to  obtain  spectral  information  at  higher  frequencies,  the 
spectra  of  the  axial  and  radial  velocities  during  a  single  time  period 
(Tsl.34  sec)  were  calculated  and  averaged  over  100  cycles.  The  spectra 
at  two  radial  positions  are  shown  in  fig. 3-23:  one  on  the  oenter  line 

B  K 

(Fig. 3-23  )  and  one  taken  at  r/R*0.?3  (Fig. 3-23  ).  The  lowest  frequency 
of  the  spectra  corresponds  to  the  frequency  of  the  pulsations,  i.e. 
f=0.74  Hz. 

The  value  of  the  leading  ooeffioient  of  the  axial  velocity  compo¬ 
nent  power  spectrum  which  corresponds  to  the  forcing  frequency  grows 
witn  increasing  amplitude  of  forcing,  but  the  rest  of  the  speotrum  is 
not  influenced  by  pulsations.  The  changes  at  the  frequenoy  of  pulsa¬ 
tions  are  assumed  to  result  from  the  experimental  error  discussed  ear¬ 
lier  in  the  case  of  long  time  spectra.  Turbulent  energy  is  increasing 
towards  the  wall,  as  expected;  most  of  the  turbulent  energy  is  conta¬ 
ined  in  frequencies  lower  than  100  Hz;  the  value  of  power  speotrum 
coefficients  at  frequencies  higher,  than  500  Hz,  is  at  least  four  orders 
of  magnitude  below  the  ooeffioients  corresponding  to  the  energy  contain- 
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ing  eddies,  so  tnat  the  sampling  frequency  of  about  l.b  khz,  used  in  the 
experiment,  is  adequate. 

ins  spectra  in  Pig .3.23  are  plotted  In  log-log  coordinates.  inis 
presentation  is  necessary  because  of  the  wide  range  of  parameters,  but 
it  conseals  the  relative  contribution  of  different  frequencies  to  the 
total  turbulent  energy.  The  connection  between  the  spectral  coeffi¬ 
cients  c(f)  and  turbulent  energy  £  may  be  presented  in  the  following 
fora 


fi«  c(f)df*  f*c(f)d(logf)  3.2 

It  is  clear  from  this  presentation,  that  when  logf  is  used  instead  of 
frequency,  the  contribution  to  the  total  turbulent  energy  at  the  fre¬ 
quency  f  is  proportional  to  f*o(f).  The  Fig. 3.23°  gives  the  spectrum  of 
the  turbulent  energy  of  the  longitudinal  velocity  fluctuations  on  the 
center  line  of  the  pipe,  as  given  in  Flg.3.23a,  but  instead  of  logc(f) 
the  logarithm  of  the  produot  f*o(f)  is  given  on  the  ordinate.  Fig. 3. 23d 
presents  the  spectra  of  radial  velooity  fluctuations,  calculated  in  the 
sane  way. 

It  can  be  seen  from  the  Figure,  that  the  main  contribution  to  the 
total  turbulent  energy  of  axial  velocity  is  at  frequencies  corresponding 
to  a  btrounal  number  baaed  on  mean  bulk  velooity  IT  and  diameter  of  the 
pipe  D  (Staffl/U-)  whioh  is  of  the  order  unity  (Stal  at  fa 55  Hz).  The 
frequency  of  energy  containing  eddies  lnoreases  near  the  oenter  of  the 
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pipe.  lne  characteristic  frequency  cf  t.u  rnergy  containing  eddies  in 
tne  raui.il  iii  action  is  about  lw  -nr,  and  seems  tc  oc  less  dependent  on 
rcuiti  position. 

2 

witn  po*er  spectrum  c(f)  Known,  dissipation  spectrum  f  c(f)  and 
(  2 

dissipation  integral  jf  c(f)df  can  oe  calculated,  'me  log-log  plot  of 
dissipation  spectra  on  the  center  line  of  tne  pipe  for  unforced  flow  and 
for  two  amplitudes  of  forced  pulsations  at  T=1.;4  sec  (U{/i7s20i»  and 
Uj/7*35>)  is  presented  in  fig.3<25.  The  maximum  dissipation  region  is 
approximately  at  f=150  Hz.  The  difference  between  the  regions  of  energy 
containing  eddies  and  maximum  dissipation  is  not  very  large,  as  can  be 
expected  for  low  iteynolds  numbers  considered.  Tne  dissipation  spectra 
are  practically  uneffected  oy  pulsations,  and  tne  slight  increase  in  tee 
value  of  the  dissipation  integral  with  Increasing  results  from  changes 
in  spectra  at  the  nigh  frequency  (f>200  Hz). 

The  spectra  on  figures  3.23  give  information  which  is  averaged  over 
the  wnole  period,  it  was  mentioned  earlier,  nowever,  that  in  the  pul¬ 
sating  flow  the  turbulent  activity  is  phase  dependent.  It  was  decided 
therefore  to  divide  each  period  into  eight  equal  parts,  and  calculate 
the  spectrum  for  each  part  separately,  bach  record  thus  included  25o 
measured  data  points  and  thus  no  spectral  information  could  be  ootalned 
about  frequencies  lower  than  tt/T«6  Hz.  The  spectral  resolution,  which 
depends  on  the  duration  of  the  total  record  is  also  reduced  by  the  suo- 
dlvlsion  of  the  period,  figures  3.24a  and  3.24b  show  three  out  of  eight 
possible  speotral  distributions  at  the  initial  pressure  pnases  of  *0°, 


p-«*e  bo 


1o0°  and  270°  ^iu  for  two  amplitudes  of  bulx  velocity  pulsations  (13o» 
on  r'io0.24‘~  uno  oo#  on  Fig. 3. 24°).  The  phase  angles  chosen  snow  tne 
maximum  deviation  l y0°  and  2?0°)  of  tne  short  duration  spectra  from  long 
time  average::  spectrum  wnicn  also  resembles  the  spectrum  at  pnase  angle 
1o0°.  At  large  amplitudes  of  pulsations,  the  difference  between  the 
specrai  coefficients  measured  at  and  ♦=270°  for  a  given  frequen¬ 
cy  may  be  as  large  as  decade.  It  is  inter sting  to  note  that  effect  of 
pulsations  on  the  spectra  is  most  pronounced  at  high  frequencies  (20 
Hz<f< 100  Hz)  rather  than  in  the  immediate  neighbourhood  of  the  forcing 
frequency  (l.e.  at  f ~ 1  Hz). 

To  evaluate  the  influenoe  of  forcing  on  the  energy  containing  ed¬ 
dies  same  spectra  were  plotted  once  again  in  the  coordinates  log(f"c(f)) 

C  d 

vs.  log(f)  (Figures  3.24  and  3.24  ).  It  transpires  tnat  in  addition 
to  the  effect  on  the  amplitude  of  the  turbulent  fluctuations,  the  intro¬ 
duction  of  forcing  also  effeots  the  energy  containing  eddies,  as  it  was 
noticed  from  examining  a  single  event  in  Section  3>3«  At  higher  instan- 
teneous  Reynolds  numbers  (with  a  certain  phase  lag)  the  entire  spectrum 
seems  to  shift  towards  higher  frequencies  which  are  appropriate  pernaps 
to  the  instantaneous  Reynolds  number. 

This  behaviour  is  even  more  pronounced  on  tne  short  duration  dissi¬ 
pation  spectra,  calculated  from  the  power  speotra  and  presented  in 
Fig. 3 . 2b.  The  values  of  the  dissipation  spectral  coefficients  differ  by 
two  orders  of  magnitude  at  different  initial  phase  angles.  The  region 
of  maximum  dissipation  also  depends  strongly  on  the  phase  angle. 
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3.3.4  Dissipation  time  scale 


Tne  multiple  spectra,  presented  in  the  previous  Section,  show  that 
the  frequency  of  turbulent  fluctuations  is  phase  dependent.  They  can¬ 
not,  however,  provide  "instantaneous"  information  about  frequency,  but 
rather  an  average  over  1/6  of  the  period.  An  alternative  way  to  obtain 
phase  looked  information  on  the  characteristic  frequency  of  the  turbu¬ 
lent  fluctuations  is  to  oaloulate  the  ensamble  averaged  values  of  the 
Eulerlan  dissipation  time  scale  T£,  defined  as 


fg*1/Tg=»Ou,/3t)^/2u,‘ 


(tiinze,  p . 45 ) .  Although  Townsend  (1956)  was  very  sceptical  about  the 
phisical  meaning  of  the  so  called  dissipation  parameters,  they  still  may 
serve  as  an  indicator  defining  the  most  rapid  turbulent  fluctuations  in 
the  flow. 


The  radial  distribution  of  time  averaged  values  of  fg  for  three  am¬ 
plitudes  of  pulsations  at  T«1.34  sec  is  presented  in  Fig. 3.27*.  No  sig¬ 
nificant  influence  of  pulsations  oan  be  noticed  up  to  the  maximum  ampli¬ 
tude  used  (i.e.  U^/U  30*).  The  average  fg  attains  a  maximum  on  the 
center  line  of  the  pipe.  The  obtained  maximum  value  of  fg  (about  400 
Hz)  supplies  one  more  indication,  that  the  sampling  frequency  ohoaen  in 
this  investigation  (about  1500  Hz)  is  quite  adequate.  The  value  of  time 
averaged  fg  increases  sharply  with  deareasing  distance  from  the  wall. 
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The  radial  distribution  of  the  amplitude  of  the  oscillations  in 
<fg>  (see  Fig.3.27*5)  resembles  the  distribution  of  the  amplitude  of 
phase  averaged  velocity  oscillations.  In  the  central  region  of  the  pipe 
the  radial  distribution  of  fg  is  approximately  uniform,  as  is  the  case 
for  the  amplitude  of  <u>,  the  difference  occurs  in  the  vicinity  of  wall, 
where  the  reduction  in  the  amplitude  of  fg  is  steeper.  The  amplitude  of 
oscillations  in  <f£>  is  roughly  proportional  to  that  of  <u>. 

The  radial  distribution  of  the  phase  lag  of  <f£>  behind  the  pres¬ 
sure  (Fig. 3.27°)  resembles  that  of  $u(r),  with  the  phase  lag  increasing 
towards  the  wall.  Die  phase  differenoe  between  <u>  and  <fg>  is  about 
20°  at  all  radial  positions,  and  is  praetioally  independent  of  ampli¬ 


tude. 
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CHAPTER  4 


DISCUSSION  AND  THEORETICAL  CONSIDERATIONS 


4.1  The  balance  of  forces  in  a  pulsating  flow 


In  a  steady  fully  developed  pipe  flow  there  are  two  types  of 
forces,  acting  on  every  element  of  fluid  and  balancing  one  another: 
pressure  forces,  resulting  from  favorable  pressure  gradient  in  the  di¬ 
rection  of  the  flow,  and  shear  forces,  which  are  caused  by  friction  on 
the  walls  and  oppose  the  motion.  In  the  case  of  non-steady  flow  a  third 
force  is  added:  an  inertia  force.  At  any  instant,  all  three  forces 
have  to  be  balanced,  forming  triangle  of  forces  for  each  frequency  in 
the  Fourier  expansion.  The  leading  term  in  the  expansion  is  considered 
and  the  presentation  of  (1.3)  is  used  (e.g.  the  analysis  refers  to  the 
frequency  of  the  pulsations  only). 


Eq.(I.IO)  may  be  rewritten  in  the  following  way: 


3<u> 

3t 


I*E>  .  i*  (r<T(r)>) 

p  3x  pr  3r  v  v  ;  J 


4.1 


where  <t>  is  the  ensemble  averaged  pulsating  part  of  the  shear  stress 


t 
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<T> 


u 


3<u> 

TT  ‘ 


p<u'v'> 


4.2 


Integrating  Eq.(4.1)  from  the  center-line  to  another  radial  position  r 
yields  after  some  manipulation 


7  It  fpr'<u(r,)>dr'  =  -  +  ~~r^—  4.3 
o' 

The  Integral  on  the  left  hand  side  represents  an  instantaneous  rate  of 
mass  flow  through  the  central  region  of  the  pipe  up  to  the  radial  posi¬ 
tion  r,  and  will  be  designated  further  p Q(r) .  Making  use  of  Eq.(1.3) 
and  considering  the  fundamental  frequency  only,  the  differential  equa¬ 
tion  (4.3)  reduces  to  an  algebraic  expression 


-i<J>  (r) 

—  pQ,(r)e  * 

1 


1  3px  (r)  i4»x(r) 

2  3x  r  6 


4.4 


bq. (4.4) ,  being  a  complex  one,  represents  in  fact  two  independent  equa¬ 
tions,  for  the  real  and  for  the  imaginary  part.  When  tne  experimental 
information  about  flow  velocities  and  pressures  is  available,  Eq.(4.4) 
contains  two  unknown  values:  the  amplitude  of  the  shear  stress  t  ^  and 
its  phase  angle  $x.  Both  values  can  thus  be  determined  from  equations 
(4.4). 


An  alternative  approach  to  Eq.(4.4)  is  to  regard  it  as  a  veotor 
equation  representing  a  balance  among  three  vectors:  QsQ1(r)exp(-i^^) , 
?»0.53P,/3x  and  x  «(x1/r)exp(i$T) ,  where  ♦  and  $T  represent  the  phase 
angles  of  mass-flux  and  shear  relative  to  the  pressure.  In  turbulent 
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flow  the  mass  flux-lags  behind  the  pressure  by  an  angle  ranging  from  0° 
to  90°,  depending  on  frequency  and  mean  Reynolds  number  (see  Sections 
3.2.2  and  3.2.4).  The  fact  that  the  mass  flux  term  on  the  left  hand 
side  of  Eq. (4.4)  is  multiplied  by  i  means  that  an  angle  of  90°  is  added 
to  this  phase  angle  in  order  that  the  mass  flux  term  on  the  left  hand 
3ide  of  Eq.(4.4)  will  lead  the  pressure  term.  The  following  sketch 
shows  qualitatively  the  relevant  vectors  and  phase  angles. 


It  is  seen  from  the  sketch,  that  both  amplitude  ^  and  phase  angle  <|>T  of 
the  shear  stress  vector  are  critically  dependent  on  the  relative  length 
of  P  and  Q  and  on  the  phase  angle  between  them,  consequently  even  a  re¬ 
latively  small  error  in  one  or  both  quantities  which  are  measured  exper¬ 
imentally,  can  lead  to  a  significant  error  in  the  shear  stress  vector. 

In  laminar  pulsating  flow  the  radial  distribution  of  the  shear 
stress  has  been  calculated  theoretically  (Uchida  (1956)),  and  the  solu¬ 
tion  for  the  amplitude  of  <t  >,  normalized  by  the  length  of  the  pressure 
vector,  is  represented  by  a  solid  line  in  Pig. 4.1  for  a»9.  (Tsl.3 
sec).  The  slope  of  the  curve  in  the  central  region  of  the  pipe  is  very 
small,  and  it  differs  notably  from  zero  only  inside  the  Stokes  layer. 
This  behaviour  of  the  radial  dependence  of  t  ?/r  in  the  laminar  flow 
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follows  from  the  fact,  that  in  the  central  region  of  the  pipe  the  accel¬ 
eration  term  is  practically  in  phase  with  pressure,  and  both  terms  oal- 
ance  one  another.  The  viscous  term  Decodes  more  important  in  the  stokes 
layer,  wnere  tne  radial  gradients  in  the  amplitude  and  pnase  of  velocity 
pulsations  increase  (rig.j.d  and  rig. 3.5).  The  increased  amplitude  of 
tne  velocity  oscillations  renders  a  longer  acceleration  vector;  the 
changes  in  the  veLocity  phase-lag  increase  the  angle  between  the  accel¬ 
eration  vector  and  the  pressure  vector;  thus  the  shear  stress  term 
which  closes  the  force  triangle  increases,  near  the  wall  the  amplitude 
of  velocity  oscillations  decreases  causing  the  acceleration  vector  to 


become  shorter  than  the  pressure,  and  the  phase  lag  of  the  shear  stress 
behind  the  pressure  tends  approximately  to  45°.  The  following  sketch 
shows  qualitatively  (but  not  to  scale)  the  force  triangles  in  laminar 
flow  at  three  radial  locations:  in  the  central  region,  in  the  Stokes 


center(r/R=0.24) 


Stokes  layer 
(r/R=0.85) 


wall 

(r/R=l) 


The  relative  length  of  the  shear  stress  vector  in  laminar  flow  is 


small  with  the  exception  of  the  wall  region  and  consequently  the  pres¬ 
sure  and  acceleration  are  collinear  and  equal .  This  identity  provides 
an  opportunity  to  check  the  consistency  of  the  measurments  and  oorrect 
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the  pressure  vector  which  in  all  probability  contains  the  largest  error. 

-*•  -»  2 

The  measured  values  of  P  and  i  in  tne  central  region  of  the  pipe 
in  laminar  flow  were  equated,  and  a  coefficient  of  the  correction  for 
the  length  of  tne  pressure  vector  was  determined.  This  coefficient  was 
used  to  correct  the  length  of  the  pressure  vector  in  the  turbulent  flow 
at  otherwise  identical  parameters  (i.e.  fie,  T,  piston  displacement). 
The  value  of  the  corection  coefficient  was  usually  about  1.15.  It  is 
worth  noting  here,  that  the  shear  stress  term  appears  in  the  form  T/r. 
This  term  does  not  vanish  in  the  center  of  the  pipe,  contrary  to  the 
shear  stress  itself,  and  therefore  there  is  no  way  to  correct  the  rela¬ 
tive  length  of  the  other  two  experimental  vectors  using  the  known  value 
of  TsO  in  the  center  of  the  pipe. 

The  shear  stresses  calculated  from  Eq.(4.4)  with  the  aformentioned 
correction  factor  for  the  pressure  can  be  checked  experimentally  by  com¬ 
parison  with  the  measured  Reynolds  stresses.  Pig. 4.1  represents  the 
calculated  radial  distribution  of  T^/p*  normalized  by  the  modulus  of  the 
pressure  vector  P,  as  compared  with  the  theoretically  calculated  laminar 
distribution.  Pour  different  amplitudes  of  pulsations  are  shown  in 
Pig. 4.1  for  the  period  T=1.34  sec. 

An  identical  way  of  normalization  was  used  in  Section  3-3.2  for  the 
radial  distribution  of  the  amplitudes  of  pulsations  of  Reynolds  stresses 
shown  in  Fig. 3. 20.  The  vector  diagrams  provide  a  reason  for  it.  As  it 
was  mentioned  in  Section  3*2,  the  amplitude  of  velocity  pulsations  is 
proportional  to  the  amplitude  of  pressure  for  both  laminar  and  turbulent 
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flows,  as  long  as  this  amplitude  is  not  too  high,  $u  is  practically  in¬ 
dependent  of  amplitude.  This  implies,  that  the  force  triangles  for  dif¬ 
ferent  amplitudes  of  pulsations  at  a  constant  frequency  are  similar,  and 
thus  the  amplitude  of  pulsations  of  the  shear  stress  is  also  proportion¬ 
al  to  the  pressure. 

It  is  seen  from  Fig. 4.1,  that  this  way  of  normalization  leads  to 
radial  distribution  of  which  collapse  fairly  well  on  the  same  curve 
provided  the  very  large  amplitudes  of  forcing  are  excluded  (U./U&35}). 
The  comparison  of  the  measured  Reynolds  stress  with  callculated  ac¬ 
cording  to  Eq.(4.4)  is  fairly  satisfactory.  In  the  central  region  of 
the  pipe  9<u>/  3r=0  (see  Section  3.2.1),  thus  <t  >/p  *-><u'v,>,  yet  the 
measured  <u'v,>1  are  consistently  below  T ^ .  The  maximum  discrepancy 
however,  is  less  than  20?  (Fig. 4.1)  thus  giving  an  estimate  of  the  error 
in  the  use  of  Eq.(4.4).  If  one  does  not  correct  the  pressure  vector  the 
discrepancy  between  the  measured  and  calculated  t1  incresases  slightly 
without  affecting  the  qualitative  results. 

Fig. 4. 2  presents  the  radial  distribution  of  the  amplitude  of  the 
pulsations  of  the  shear  stress,  normalized  in  the  same  way  as  in  the 
previous  Figure,  for  four  different  frequencies.  The  relative  lengths 
of  the  shear  stress  vector  and  of  the  pressure  vector  are  obviously 
strongly  dependent  on  frequency.  At  low  frequency  the  length  of  T  at 
the  wall  nearly  equals  to  that  of  I,  as  it  is  in  steady  flow,  and  the 
contribution  of  the  acceleration  term  is  manifested  mainly  in  the  phase 
shift  between  t  and  It  is  natural  to  assume,  that  at  frequencies 


Page  75 


lower  thar<  the  ones  measured  this  phase  lac  eventually  vanisnes,  and  tne 
behaviour  of  tne  pulsating  turbulent  pipe  flow  aoes  not  differ  from  tne 
steady  flow,  ine  length  of  x  relative  to  the  amplitude  of  the  pressure 
pulsations  decreases  with  increasing  frequency,  and  the  behaviour  of  tne 
flow  in  the  central  region  of  the  pipe  resembles  the  laminar  case,  where 
the  influence  of  the  shear  stress  is  negligible. 

It  should  be  noted  however  that  the  amplitude  of  the  shear  stress 
does  not  change  notably  with  frequency  wnen  normalized  by  the  flow  rate 
amplitude  Qj  rather  than  the  pressure.  This  is  also  the  case  for  the 
RMS  values  of  the  longitudinal  velocity  oscillations  (see  Pig. 3. 17).  It 
may  be  concluded  that  the  oscillations  in  shear  stress  and  RMS  values  of 
velocity  fluctuations  are  functions  of  the  oscillating  Reynolds  number. 

There  exists  a  phase  difference  between  the  oscillations  in  mass 
flux  and  in  shear  stress.  The  phase  lag  of  <x  > 4>r  as  calculated  from 
the  £q.(4.4),  shows  only  slight  dependence  on  amplitude  (Pig. 4. 3),  with 
the  possible  exception  at  very  high  amplitude  of  forcing  (11^/1^35^).  In 
the  latter  case  increases  notably.  The  dependence  of  the  phase  angle 
on  frequency  is  considerable  (Fig.  4.4).  The  phase  angle  <t>T  calculated 
from  Eq. (4.4)  and  the  angle  $uy  which  was  deduced  from  the  measured  Re¬ 
ynolds  stresses  are  in  very  good  agreement  (Fig. 4.3).  The  difference 
between  the  two  quantities  is  less  than  10°.  The  phase  angle  between 
the  pressure  and  the  shear  stress  increases  with  increasing  frequency. 
♦T  decreases  near  the  wall,  this  deorease  is  ouch  more  pronounced  at 


higher  frequencies 
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The  following  sketch  (which  can  not  be  drawn  to  scale)  shows  rela¬ 
tive  lengtn  and  directions  of  the  relevant  vectors  at  a  given  radial  lo¬ 
cation  (r/R=0.o1)  for  tnree  frequencies  (all  vector  moduli  are  normal¬ 
ized  by  the  modulus  of  P),  as  calculated  from  experimental  data  in  ac¬ 
cordance  with  l£q.(4.4): 


T  T 


Increasing  the  frequency  results  in  smaller  phase  angles  between  the 
pressure  and  the  acceleration  vector.  The  length  of  Q  relative  to  the 
pressure  decreases  with  increasing  frequency  (in  analogy  to  the  behavi¬ 
our  of  cj^  a)  in  laminar  flow  (Fig. 2. 2)),  but  the  product  ioQ1  which 
corresponds  the  the  acceleration  vector  increases  with  increasing  fre¬ 
quency.  Thus,  while  examining  the  instantaneous  force  balance,  t  oeing 
related  to  the  acceleration  should  be  compared  with  a  characteristic  oft 
term,  while  in  examining  the  influence  of  forcing  on  the  turbulent  in¬ 
tensities  or  Reynolds  stress  should  be  compared  with  a,  characteristic 

•¥ 

amplitude  of  velooity  and  henoe  with  Q^.  The  phase  lag  of  T  behind  P 
increases  with  frequency,  as  shown  in  Fig. 4. 4. 

At  a  given  radial  position  there  is  no  practical  difference  between 
<t  >/pand  -<u'v*>.  The  phase  angle  of  <  t>  may  thus  be  considered  to 
be  equal  to  the  phase  angle  of  the  Reynolds  stresses.  The  phase  lag  of 
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tne  Reynolds  stre33es  behind  tne  mass  flux  tnus  ranges  from  nearly  zero 
at  low  frequencies  to  approximately  o0°  at  T=0.7o  sec.  * t  nigh  frequen¬ 
cies  Reynolds  stress  is  nearly  in  antiphase  with  tne  acceleration  vec¬ 
tor.  since  the  radial  dependence  of  the  phase  angle  of  the  longitudinal 
velocity  fluctuations  is  similar  to  the  dependence  of  the  Reynolds 
stresses  (because  of  production  term),  it  follows  from  the  force  trian¬ 
gles,  that  at  high  frequencies  the  acceleration  is  associated  with  re¬ 
duction  in  the  turbulent  activity.  This  fact  is  widely  recognized  in 
the  literature  (e.g.  Narasiaha  and  Sreenivasan  (1979))  and  a  causality 
of  relaminarisation  by  acceleration  was  suggested,  but  the  present  ex¬ 
perimental  results  suggest  a  different  interpretation. 

At  long  periods  of  pulsations  (Ts4.4b  sec),  the  rate  of  change  of 
the  oulk  velocity  is  small  enough  for  the  turbulent  structure  to  accomo¬ 
date  itself  to  the  instantaneous  value  of  the  phase-loctced  average  velo¬ 
city  distribution.  The  changes  in  the  turbulent  characteristics  of  the 
flow  sure  therefore  in  phase  with  the  mass  flux  and  with  the  pulsating 
Reynolds  number.  The  situation  resembles  laminar  pulsating  pipe  flow, 
where  at  a<1,  the  velocity  at  any  instant  has  a  Poiseuille  distribution 
corresponding  to  the  instantaneous  value  of  the  pressure  gradient.  The 
period  of  pulsations  for  whioh  the  turbulent  structure  responds  to  the 
instanteneous  mean  flow  (T»4. 5  sec)  is  still  very  short  in  terms  of  the 
frequenoy  parameter  (  a«4.9  rather  than  1).  Thus  the  characteristic 
responoe  time  of  turbulent  fluctuations  to  a  change  in  mean  flow,  is 
much  shorter  than  the  time  needed  for  the  laminar  flow  to  acoomodate  it¬ 
self  to  the  ohanging  pressure  gradient. 
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when  the  frequency  of  pulsations  increases,  the  phase  lag  of  the 
mass  flux  relative  to  the  pressure  tends  to  90°,  and  all  tnree  vectors 
in  the  force  triangle  tend  to  be  coll inear.  This  implies  that  the  Rey¬ 
nolds  stresses  lag  behind  the  mass  flux  by  90°  and  cannont  exceed  this 
angle  under  any  condition.  At  high  enough  frequencies  the  turbulent  ac¬ 
tivity  is  minimum  when  the  acceleration  is  maximum  because  the  relative 
contribution  of  the  shear  stresses  to  the  total  balance  of  forces  decre¬ 
ases.  The  "laminar! zing"  effect  of  the  acceleration  is  thus  obtained  in 
the  limiting  case  of  very  rapid  oscillations. 


4.2  The  Eddy  Viscosity  Model  in  the  Time  Dependent  Flow 


In  pulsatirg  pipe  flow  amplitudes  and  phase  angles  of  the  oscillat¬ 
ing  pressure,  velocity  and  shear  stress  are  mutually  dependent. 
Therefore  the  qualitative  difference  in  the  radial  distribution  of 
in  laminar  and  turbulent  flows  (Fig. 3. 5)  results  from  the  different  be¬ 
haviour  of  the  oscillating  part  of  the  shear  stress  in  both  flows.  The 
oscillating  parts  of  the  laminar  and  turbulent  velocity  profiles  at  high 
values  ofxare  vaguely  similar,  suggesting  that  a  proper  use  of  an  eddy 
viscosity  model  may  provide  a  qualitative  answer  for  the  behaviour  of 
the  oscillating  turbulent  flow,  including  in  particular  the  vastly  dif¬ 
ferent  radial  distrubution  of  the  phase.  The  eddy  visoosity  model  has 
to  be  modified  in  order  to  provide  a  qualitatively  correct  description 
of  the  time  dependent  turbulent  flow. 


Page  79 


The  suggested  expression  of  the  kinematic  eddy  viscosity  for  the 
time  dependent  portion  of  the  flow  is  based  on  the  steady  flow  model  of 
Van-Driest  (195o)  with  the  correction  of  the  wake  function  of  Coles 
(1954),  as  proposed  by  Yahalom-Dioant  (1974).  The  kinematic  eddy  vis¬ 
cosity  eQ  for  steady  flow  is  given  by 


■  m 


4.5 


The  velocity  gradient  3u/3r  is  composed  of  two  parts,  the  first  repre¬ 
sented  by  Van-Driest  model : 


3ux  (r/R)u2/v 

3r  ,  R.r  ,  r  (R-r)u*  7~ 

l+[1+4k2(ii)u^  £  (l-exp( - —)))’ 

VA 

for  which  the  empirical  constants  are:  A+=26;  k=0.4;  and  the 

part  is  the  wake  correction  for  the  central  region  of  the  pipe: 


4.6 


second 


3u„ 


[l-exp(-2Ru,/vA*)sin(ir(l-  £)] 


4.7 


where  the  additional  empirical  constant  II  was  chosen  to  be  11x1.02. 


The  eddy  viscosity  for  the  time  dependent  part  of  the  flow  is  dif¬ 
ferent  from  that  for  the  steady  flow.  Following  Maruyama  (1974)  and 
using  the  Prandtl's  mixing  length  approach,  the  expressions  for  the 
shear  stress  and  the  eddy  viscosity  are: 


P£, 


3u 

17 


e 


o 


*  pi.2 


3u 

17 


4.6 
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where  l  is  tne  mixing  lengtn.  introducing  tne  decomposition  for  tne 
steady  and  time  dependent  parts  of  the  pnase  averaged  velocity  u=u+<u>, 
suostituting  it  into  £q.(4.o),  and  neglecting  terms  of  seoond  order 
yields: 

'  ■  ♦  *»*  §f  ^  M 

thus  implying  for  the  eddy  viscosity  of  the  oscillating  flow  e, 
el  =  2pZl  fr  =  2eo  3 (W8r)  4'9i 


The  simple  model  of  Prandtl  which  uses  a  constant  length  scale,  gives 
for  the  time  dependent  part  of  the  shear  stress  an  eddy  exchange  coeffi¬ 
cient  ^  which  is  twice  the  conventional,  steady  eddy  viscosity. 

For  a  more  complex  model,  like  the  one  given  in  equations  (4.5) 
through  (4.7),  the  relation  between  eQ  and  e 1  is  not  as  simple.  In 
order  to  obtain  the  eddy  viscosity  which  represents  the  dependence  of 
the  shear  stress  on  the  velocity  profile,  the  radial  distributions  of 
the  velocity  gradient  and  of  the  shear  stress  was  calculated  for  two 
slightly  different  values  of  the  friction  velooity  u#  from  equations 
(4.6)  and  (4.7).  The  eddy  visoosity  for  the  oscillating  flow  was  then 
calculated  following  (4.9a)  from  the  equation: 

(r/R)AuJ 

ei  -  L(WIt)  4•1, 
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Fig. 4.5  presents  so  calculated  radial  distrioution  of  tfte  oscillating 
odcy  viscosity  as  compared  to  tne  steady  one. 

It  was  tacitly  assumed  in  this  derivation,  that  the  shear  stress  is 
corresponding  to  the  instantaneous  value  of  the  local  velocity  gradient, 
irrelevant  to  the  rate  of  change  of  the  phase-averaged  flow.  The  care- 
full  examination  of  the  experimental  data  however  reveals,  that  this  is 
not  always  true.  From  the  Eq.(4.2)  and  the  definition  of  the  eddy  vis¬ 
cosity  one  obtaines 

(erv)  ^5? a  -<u'v,>  4.11 


The  phase  angle  of  the  Reynolds  stresses  is  thus  related  to  the  velocity 
gradient  3<u >/3r.  With  the  radial  distribution  of  both  amplitude  u^r) 
and  phase  angle  of  velocity  oscillations  $u(r)  Known,  the  time  depen¬ 
dent  velocity  gradient  is: 


3u  (r)  3$  (r)  i(wt*A  (r)) 

-17 - -ST-'* 


4.12 


The  phase  angle  of  3<u>/3r  may  thus  be  calculated  from  Eq.(4.12).  The 
experimentally  observed  radial  dependencies  of  uf(r)  and  <t’u( r ) ,  aa  pre¬ 
sented  in  Section  3.2,  are  very  slowly  varying  functions  of  r  with  the 
exception  of  the  wall  area,  thus  the  differentiation  of  these  curves 
will  inevitably  cause  a  very  large  experimental  error.  Consequently, 
the  accurate  calculation  of  the  eddy  visoosity  for  the  oscillating  part 
of  the  flow  from  the  experimentally  measured  values  of  <t  >  and  3  <u>/3r 
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j^ons  impossiole  in  practice. 

inere  exists,  nowever,  a  point,  at  whicn  the  pnase  lag  $u(r)  at~ 
tains  a  minimum,  so  that  3$u(r)/3r=Q  at  this  radial  position,  and  the 
pnase  angles  of  3<u>/3r  and  <u>  are  identical  (£q.(4.12)).  The  experi¬ 
mental  data  of  Section  3*3.2  show  a  difference  in  phases  between  <u>  and 
<u*v*>  at  this,  as  well  as  at  other  radial  positions.  It  follows  there¬ 
fore  from  fcq.i4.11),  that  the  turbulent  part  of  the  eddy  viscosity  for 
the  time  dependent  flow  has  to  be  a  complex  function  which  introduces  a 
phase  shift  between  3<u>/3r  and  <u*v'>. 

The  complex  nature  of  the  eddy  viscosity  for  the  oscillating  flow 
can  oe  related  to  the  finite  relaxation  time  in  turbulent  flows  (Nee  and 
Kovasznay  (1909),  Narasimha  and  Praohu  (1972)).  The  concept  of  the  re¬ 
laxation  time  refers  to  the  time  period  necessary  for  the  turbulent 
structure  to  adjust  itself  to  a  step-function  change  in  the  mean  flow, 
bee  and  Kovasznay  and  Prabhu  and  Narasimha  (1972)  proposed  model  equa¬ 
tions  for  calculating  the  shear  stress  in  a  flow  with  non-stationary 
mean  values,  based  on  the  diffusion  equation.  A  different  aproach  will 
be  taxen  here.  The  equation  relating  the  shear  stress  to  velocity  gra¬ 
dient  is  modified  in  order  to  include  a  responce  time  of  the  Reynolds 
stress : 


<u ' v  *  > 


(e^v) 


3<u> 


Q  3<u'v'> 

-0— 5t — 


4. 13 


The  second  term  on  the  right  hand  side  of  the  £q.(4.13)  represents  the 
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"inertia"  of  the  Reynolds  stress,  whicn  nay  be  large  whenever  the 
cnanges  in  <u'v’>  are  rapid  in  comparison  witn  characteristic  time  9  , 
but  vanishes  for  infinitely  slow  processes.  Substituting  into  the 
i£q . ( 4 - 13)  the  presentation  (1.3)  for  the  time  dependent  part  of  the  Rey¬ 
nolds  stress  <u'  v'>=(u'  v* )  ^expdjs  t+  4>uv),  one  obtains: 


<u.v.>  =  El~V,  3^u>  =  £l'V  -iarctg(u9)  3<u> 

l+iu>0  3r  -  9r 

/i+u202 


4. 14 


which  gives  an  expression  for  the  oscillating  eddy  viscosity: 


1  .-iarctK((u0) 

'  e 

/l+(l)202 


4.15 


£q. (4. 15)  indicates  that  the  argument  of  the  complex  eddy  viscosi¬ 
ty,  as  well  as  its'  absolute  value  are  frequency  dependent.  At  high 
frequencies  relative  to  the  relaxation  time  0  (  u)0>>1)  the  absolute 
value  of  the  eddy  viscosity  for  the  oscillating  flow  decreases,  and  the 
Reynolds  stresses  become  independent  of  the  phase  angle.  Thus,  the  tur¬ 
bulence  is  "frozen"  and  depends  on  the  mean  flow  only.  The  pulsations 
are  governed  by  the  molecular  viscosity  v  only,  and  are  in  this  sense 
laminar,  while  the  mean  flow  is  fully  turbulent. 

As  mentioned  earlier,  the  phase  difference  between  <u'v'>  and 
3<u>/  3r  is  known  accurately  at  one  radial  position  only,  where  the 
phase  of  the  velocity  gradient  equals  to  the  phase  of  the  velocity  it¬ 
self.  The  experimentally  measured  $u  at  Tsi.34  sec  and  Uj/U«20}  at- 


i 


i 
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tains  aininu<;)  at  r/R=0.72  (fig. 3. 1b),  and  the  phase  difference  between 
the  velocity  and  the  Reynolds  stress  is  at  this  radial  position  30° 
(fig. 3.21).  Tne  relaxation  time  0  ,  obtained  from  the  Eq.(4.15)  is 
therefore  8  =(  ir/6)/(2  ir/T)=0.12  sec.  Although  the  relaxation  time 
should  in  general  depend  on  the  radial  position,  it  is  convenient  to  as¬ 
sume  it  to  be  constant  across  the  pipe.  The  relaxation  time  therefore 
is  a  universal  constant  for  a  given  mean  Reynolds  number.  Narasimha  and 
Prahbu  (1972)  used  a  similar  approach  to  calculate  the  influence  of  the 
finite  relaxation  time  on  the  development  of  the  turbulent  wake  in  res- 
ponce  to  a  steep  change  in  pressure  gradient.  Neglecting  the  dependence 
of  the  relaxation  time  on  the  transverse  coordinate  proved  in  their  case 
to  have  no  significant  effect  on  the  results  of  the  calculations  nor  on 
the  agreement  with  experimental  results. 

Following  equation  (4.15),  one  may  define  a  critical  frequency  of 

pulsations  f  /2iTs1/2jr6  si. 3  Hz.  At  w  <  u)  the  modulus  of  g  ,  is 

cr  cr  cr  i 

practically  independent  of  frequency  and  equals  to  the  value  of  an  eddy 

exchange  coefficient  for  very  slow  pilsations.  At  the  oscillat¬ 

or 

ing  eddy  viscosity  decreases  to  1/^2  of  its  value  at  (#+  0.  When  u>  >1/0, 
2 

and  (u8  )  »1,  lej'l/w  ,  and  thus  the  amplitude  of  oscillations  in 

shear  stress  decreases  rapidly  with  increasing  frequency  for  a  given 

value  of  3<u>/  dr.  The  corresponding  critical  period  of  the  pulsations 

is  in  this  ease  T  =1/f  a0.75  sec.  For  the  turbulent  structure  to  be¬ 
er  or 

come  frozen,  i.e.  for  the  oscillating  eddy  viscosity  to  vanish,  the 

frequency  of  pulsations  has  to  be  much  higher  than  f  ,  as  may  be  de¬ 
er 

duced  from  Eq.(4.l5).  On  the  other  hand,  even  at  f<f  the  oscillating 

cr 
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eddy  viscosity  still  introduces  a  significant  pnase  angle  between  <u'v'> 
and  9<u>/  3r,  thus  influencing  on  the  total  force  oalance  in  the  oscil¬ 
lating  part  of  the  flow. 

Ramaparian  and  Tu  (I960)  used  a  criterion,  proposed  by  Rao  et  al 
( 1 97 1 ) ;  Laufer  and  Badri  harayanan  (197D  for  tne  cnaracteristic  res- 
ponce  time  in  a  turbulent  boundary  layer  U0  /  6=5,  in  order  to  distin¬ 
guish  slow  pulsations  from  rapid  ones.  The  adaptation  of  tne  criterion 
to  the  pipe  flow  required  that  the  boundary  layer  thicicness  is  replaced 
by  a  diameter  D  of  the  pipe.  The  direct  adaptation  of  results  obtained 
in  the  boundary  layer  to  a  pipe  flow  seems  questionable,  and  the  radius 
is  probably  a  better  substitute  for  6,  than  D.  Still,  the  use  of  the 
aoovementioned  criterion  for  the  flow  parameters,  used  in  the  present 
invesigation,  gives  a  characteristic  time  0=5D/u=O.Oo  sec,  which  is  in 
fair  agreement  whitnin  an  order  of  magnitude  with  the  estimated  value  of 
0=0.12  sec  mentioned  earlier.  Ramaparian  and  Tu  concluded  that  the 

turbulent  structure  is  independent  of  the  phase  angle  at  f~  f  ,  yet  in 

cr 

the  present  investigation  experimental  results  and  model  considerations 
indicate,  that  much  higher  frequencies  are  necessary  in  order  to 
"freeze"  the  turbulence.  The  characteristic  timed,  when  calculated 
from  the  criterion  0=5D/U,  decreases  with  increasing  Re,  thus  the  phase 
shift  between  the  velocity  and  the  turbulent  structure  has  to  decrease 
with  increasing  time-mean  bulk  velocity.  The  Reynolds  number  dependence 
of  the  phase  shift  between  <u>  and  J< u'^>  (Fig.3. 15b)  confirms  this  con¬ 
clusion. 
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It  is  intersting  to  note  that  the  empirical  relation  of  Mizushina 

et  al  (1973  )  (£q.(1.12)  leads  in  our  case  to  the  same  value  of  T  =0.75 

cr 

sec.  These  authors  also  referred  to  the  paper  of  Rao  et  al  (1971); 
their  conclusion  is  however  opposite  to  the  results  ootained  in  the  pre¬ 
sent  work:  Mizushina  et  al  found  that  for  T>T  the  pulsations  do  not 

CP 

effect  the  turbulent  structure,  and  for  very  rapid  pulsations  (T<T  ) 

cr 

turoulent  intensities  become  phase  dependent.  This  conclusion  contrad¬ 
icts  to  the  physical  sense  and  experimental  results  and  seems  to  be 
error . 


The  eddy  viscosity  calculated  from  the  Eqs.(4.5)  to  (4.7),  (4.10) 
and  (4.15),  was  substituted  into  the  Navier-Stokes  equation  for  the  os¬ 
cillatory  part  of  the  flow  to  give: 


iw<u>  = 


I  5<£>  +  II 

p  3x  r  3r 


(r^O) 


3<u>  , 
3r  ' 


4.  16 


with  boundary  conditions: 


3<u> 

3r  r=0 


=  0 


=  0 


Eq. (4. 16)  was  solved  numerically  using  an  implicit  finite  difference 
scheme  by  the  Gaussian  elimination  method  (Chow  (1979),  P«74). 


The  results  of  the  calculations  are  presented  in  Fig. 4. 6,  showing 
the  radial  distribution  of  the  amplitudes  and  phase  angles  of  the  oscil¬ 
lations  in  velocity  and  shear  stress.  The  calculations  were  performed 
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for  jo an  he=4CuO  and  1=1.0  3ec.  a  qualitative  agreement  witn  tne  exper¬ 
imental  results  was  acnieved  in  tnis  way.  Tne  calculated  pna3e  lag  of 
tne  velocity  oehiria  tne  pressure  behaves  as  tne  experimentally  ooserved. 
Tne  calculated  angle  sligntly  decreases  with  the  increasing  distance 
from  the  center  line,  then  increases  near  the  wall.  It  should  be  noted 
that  tne  efforts  to  use  the  real  value  of  the  oscillating  eddy  viscosity 
resulted  in  the  calculated  radial  distribution  of  <J>  which  resembled 
<}>u  in  laminar  flow,  and  decreased  monotonically  with  increasing  r/fi. 

In  the  immediate  proximity  of  the  wall  the  rate  of  change  of 
<J>u  with  radius  decreases  sharply.  The  comparison  of  the  <f>^  which  was 
calculated  from  the  experimental  data  according  to  Eq.(4.4)  and  shown  in 
Fig. 4. 4,  with  at  the  same  large  radial  location  (snown  in  Fig. 3. 5) 

reveals,  that  the  phase  difference  between  them  is  small.  The  turbulent 
(complex)  part  of  the  eddy  viscosity  Deing  negligiole  in  this  region, 
the  assumption  can  be  made  that<T>  has  to  be  in  phase  with  3<u>/  gr, 
tnua  the  phase  angle  gradient  3$y(r)/  3r  has  to  be  small  near  the  wall 
as  it  follows  from  Eq.(4.12).  The  behaviour  of  <{>u(r)  in  the  vicinity 
of  wall  calculated  numerically  seems  therefore  reasonable,  but  no  direct 
experimental  data  is  available  in  this  region. 

The  calculated  values  of  both  ij>u(r)  and  (r)  differ  from  the  ex¬ 
perimentally  measured  angles,  implying  that  the  theoretical  model  is  not 
sopnisticated  enough  to  provide  an  exact  balance  between  viscous,  iner¬ 
tia  and  pressure  forces,  and  perhaps  is  not  valid  for  the  low  Reynolds 
number  used. 


.  eadteSftfov  i- ' 


p*se  66 


CHAPTER  5 


CONCLUSIONS 


The  pulsating  flow  in  pipe  was  investigated  experimentally, 
rieasurments  were  carried  out  in  fully  developed  laminar  and  turbulent 
flows  consecutively  while  all  other  flow  parameters  were  kept  constant; 
enabling  a  detailed  comparison  between  the  laminar  and  turbulent  flows 
at  otherwise  Identical  conditions.  Most  measurments  were  taken  at  a 
mean  Re =4000,  for  the  periods  of  pulsations  from  0.5  sec  to  5  sec  cor¬ 
responding  to  a  dimensionless  frequency  parameter  5<a=R*/w7v<15. 

By  considering  the  fundamental  forcing  frequency  only,  the  regular 
time- dependent  component  of  the  flow  parameters  is  represented  by  the 
radial  distributions  of  amplitude  and  phase-angle  of  the  oscillations. 
It  seems,  however,  interesting  and  possible  to  expand  the  analysis  to 

include  higher  harmonics,  whioh  are  important  whenever  turbulence  is 

* 

considered.  The  restriction  of  the  analysis  to  the  fundamental  frequen¬ 
cy  made  it  possible  to  represent  the  balanoe  between  the  inertia,  pres¬ 
sure  and  shear  forces  in  pulsating  flow  by  a  vector  triangle.  The  vec¬ 
tor  approaoh  is  helpful  in  understanding  the  physical  process  at  hand. 

The  time  mean  properties  of  the  flow  were  found  to  be  praotioally 
indpendent  of  pulsations  in  both  laminar  and  turbulent  regimes.  A  small 
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increase  in  tne  turoulent  energy  at  the  nigr.  frequency  enu  of  tne  spec- 
trua  was  observed  at  relatively  nigh  amplitudes  of  forcing,  no  accumu¬ 
lation  of  the  turoulent  energy  at,  or  near  tne  frequency  of  pulsations 
was  ooserved  in  this  investigation. 

Tne  theoretical  analysis  of  the  fully  developed,  laminar,  pulsating 
pipe  flow  provides  a  good  estimate  for  the  measured  quantities.  The 
theory  describing  tne  laminar  pulsating  flow  in  the  entrance  region  of 
the  pipe  (Atabek  and  Cheng  (1961))  which  uses  a  linear  approximation  for 
tne  inertial  term  in  equation  of  motion  underestimates  tne  importance  of 
this  region  in  comparison  with  tne  experimental  results. 

The  dependence  of  the  bulk  velocity  oscillations  on  pressure  was 
found  to  oe  similar  in  laminar  and  turbulent  pulsating  flows.  The  am¬ 
plitude  of  the  oulk  velocity  pulsations,  depends  linearly  on  tne  pres¬ 
sure  in  both  flow-regimes;  so  that  equal  amplitudes  of  pressure  pulsa¬ 
tions,  result  in  nearly  identical  amplitudes  of  velocity  regardless  of 
the  flow  regime.  The  pn«se-lag  of  the  bulk  velocity  relative  to  the 
pressure  is  different  in  laminar  and  in  turbulent  flows,  the  phase  lag 
in  the  turbulent  flow  being  usually  notably  smaller  tnan  in  the  corres¬ 
ponding  laminar  flow.  The  radial  dependence  of  tne  amplitude  of  the 
velocity  oscillations  was  also  found  to  be  different  in  laminar  and  tur¬ 
bulent  flows.  In  laminar  flow,  the  maximum  amplitude  of  the  velocity 
oscillations  occurs  in  the  Stokes  layer,  resulting  in  an  ri-like  shape  of 
velocity  profile.  In  the  turbulent  flow  the  amplitude  of  velocity  os¬ 
cillations  decreases  monotonlcally  from  the  center  towards  tne  wall  pro- 
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vided  that  trie  /iccous  sublayer  is  thinner  than  the  stokes  layer.  At 
him  frequency  oi  forcing  in  tne  turoulent  flo>.  regime,  wnen  tne  stokes 
layer  is  thinner  t run  tne  viscous  sublayer,  tne  radial  distrioution  of 
tne  amplitudes  of  tne  velocity  pulsations  resmoles  tne  distribution  oo- 
served  in  the  laminar  flow.  The  radial  distrioution  of  the  pnase  angle 
♦u,  changes  significantly  when  the  flow  regime  changes  from  laminar  to 
turbulent.  In  the  laminar  flow,  <t»u  is  approximately  90°  in  the  central 
region  of  the  pipe  at  ail  frequencies  considered,  and  decreases  to  about 
45°  near  the  wall.  In  turbulent  flow,  $u  is  smaller  in  the  center  than 
in  the  corresponding  laminar  case ,  and  increases  with  increasing  fre¬ 
quency.  At  all  frequencies  considered,  the  phase  angle,  $u,  in  turbu¬ 
lent  flow  is  higher  near  the  wall  of  the  pipe  than  near  the  center,  in 

contrast  to  the  distribution  of  6  in  the  laminar  now. 

Tu 

The  turbulent  structure  in  pulsating  pipe  flow  was  found  to  be 
phase  dependent.  The  amplitude  of  the  oscillations  of  longitudinal  and 
radial  velocity  fluctuations,  as  well  as  the  Reynolds  stress,  correspond 
roughly  to  the  amplitude  of  pulsations  of  tne  bulk  velocity.  The  angle 
between  the  intensity  of  the  longitudinal  component  of  the  velocity 
fluctuations  and  the  pressure  attaines  a  minimum  in  the  region  where 
most  turbulent  energy  is  produced  and  increases  towards  the  center  line. 
The  radial  distribution  of  phase  of  the  Reynolds  stresses  resembled  that 
of  the  longitudinal  velooity  fluctuations,  but  the  minimum  is  less  pro- 
nounoed.  The  phase  lag  of  the  radial  velocity  fluctuations  is  uniform 
acoross  the  pipe,  and  approximately  equals  to  the  phase  angle  of  longi¬ 
tudinal  velooity  fluctuations  in  the  oenter  of  the  pipe.  The  phase  lag 
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of  the  turbulent  quantities  behind  the  velocity  of  forcing  increases 
with  increasing  frequency  from  zero  to  about  90°.  At  high  frequencies, 
the  turbulent  activity  reduces  to  a  minimum  whenever  the  acceleration  is 
high,  the  phase  lag  between  the  two  leads  to  the  suggestion  that  the 
turbulent  structure  needs  time  to  accomodate  itself  to  the  everchanging 
mean  flow,  thus  revealing  properties  of  "memory".  Heasurments  of  phase 
between  Reynolds  stress  and  velocity  oscillations  at  high  frequency,  en¬ 
abled  evaluation  of  a  characteristic  responce-time  of  a  turbulent  struc¬ 
ture.  This  characteristic  time  was  used  in  developing  an  eddy  viscosity 
model  for  the  oscillating  part  of  the  turbulent  pulsating  pipe  flow. 
The  model  gave  an  expression  for  a  complex  eddy  viscosity  for  the 
time-dependent  component  of  the  flow,  for  which  both  modulus  and  argu¬ 
ment  are  frequency  dependent.  The  argument  of  the  eddy  viscosity 
differs  from  zero  at  relatively  low  frequencies,  and  tends  to  90°  with 
decreasing  period  of  pulsations.  The  modulus  of  the  oscillating  eddy 
viscosity  is  independent  of  frequency  whenever  the  period  of  pulsations 
is  longer  than  the  relaxation  time  of  the  turbulence,  but  decreases  at 
higher  frequencies.  The  relaxation  time  of  the  turbulent  structure  was 
evaluated  from  a  single  point. 

The  dependence  of  the  eddy  visoosity  on  the  frequency  leads  to  the 
suggestion  that  for  a  sufficiently  high  frequencies  the  modulus  of  the 
eddy  visoosity  vanishes,  and  the  turbulent  structure  becomes  Independent 
of  forcing.  The  turbulence  therefore  becomes  "frosen",  and  the  oscil¬ 
lating  pert  of  the  flow  may  be  considered  laminar  although  the  mean  flow 
is  turbulent.  These  high  frequencies  could  not  be  achieved  in  the  pre- 
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sent  experimental  setup. 

Plow  analysis  notwithstanding,  the  responce  of  the  facility,  con¬ 
sisting  of  a  valvless  piston  pump,  a  large  settling  chamber  and  straight 
smooth  pipe,  to  a  periodic  change  in  the  volume  of  the  settling  chamber, 
was  analysed  theoretically  and  experimentally.  It  was  found  that  for  a 
given  amplitude  of  oscillations  in  volume  the  resulting  pressure  oscil¬ 
lations  depend  on  frequency.  A  maximum  responce  of  the  system  to  oscil¬ 
lations  in  volume  occurs  in  both  laminar  and  turbulent  flows  at  approxi¬ 
mately  identical  frequencies.  The  resonance  type  of  behaviour  of  the 
pressure  oscillations  is  caused  by  the  phase  relation  between  the  pulsa¬ 
tions  in  the  flow  rate  and  pressure  for  the  specific  facility.  The  in¬ 
fluence  of  the  finite  sound  velocity  on  the  linearity  of  the  pressure 
distribution  along  the  pipe  was  investigated  as  well  and  a  qualitative 
criterion  was  obtained  for  the  range  of  frequencies  at  which  a  fully  de¬ 
veloped  flow  may  be  regarded  to  be  independent  of  the  axial  coordinate. 
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